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I.  INTRODUCTION 


A  large-scale  reusable  free-field  nuclear  airblast  simulator 
based  on  fuel-air  explosives  (FAE)  has  been  investigated.  A 
fuel-air  explosive  consists  of  a  combustible  atmospheric  dispersion 
of  fuel  droplets,  fuel  vapor,  or  a  combination  of  these.  Subsequent 
explosion  of  the  fuel-air  cloud  is  induced  in  the  detonative  (as 
opposed  to  def lagrative)  mode.  This  results  in  peak  explosion 
overpressures  on  the  order  of  2  MPa  (300  psi),  developed  as  the 
detonation  front  propagates  through  the  cloud. 

The  feasibility  of  an  FAE  nuclear  airblast  simulator  is 
contingent  on  developing  a  technique  to  very  rapidly  form  an 
unconfined  fuel-air  cloud  of  the  required  size.  The  usefulness  of 
the  simulator  would  also  depend  on  the  degree  to  which  a  cloud  so 
formed  and  detonated  simulates  a  nuclear  airblast,  at  least  in  the 
sense  of  producing  correct  target  response. 

In  order  to  generate  the  highest  usable  overpressure  levels 

from  a  given  fuel-air  cloud,  hemispherical  cloud  geometry  is 

desirable.  With  other  geometries,  the  airblast  emerging  from  the 
detonated  fuel-air  cloud  would  not  be  expected  to  be  point  symmetric 
in  the  cloud  vicinity.  The  cloud  formation  scheme  that  has  been 
under  investigation  for  producing  hemispherical  clouds  is  that  of 
impulsive  liquid  fuel  injection  from  a  central  dispensing  source  on 
the  ground.  The  fuel  is  forced  at  high  pressure  through  a  large 
number  of  hydraulic  nozzles  aimed  in  directions  that  are  distributed 
uniformly  throughout  the  hemispherical  volume  of  air.  It  has  been 
estimated  that  a  1  KT  simulation  would  require  a  cloud  with  a  radius 
of  71  m  (233  ft)  and  containing  85,400  kg  (188,000  Ibm)  of  fuel 
(assuming  propylene  oxide). 

Earlier  research  at  Systems,  Science  and  Software  (S3) 

addressed  various  aspects  of  the  physical,  technological,  and 

engineering/economic  feasibility  of  an  FAE/nuclear  airblast 
simulator.  These  efforts  are  discussed  in  detail  in  References  1 


In  early  experimental  work  (Ref  1),  small-scale  fuel-air 
clouds  (1/4  ton  nuclear  equivalent,  nominal)  were  formed  and 
detonated,  and  full-scale  jets  from  single  nozzles  were  tested.  The 
small-scale  propylene  oxide/air  clouds  were  reasonably  hemispherical 
in  outline  but  they  were  not  well  filled  in  with  spray  because  too 
few  nozzles  were  used.  In  order  to  detonate  these  clouds  with  good 
results  it  was  found  necessary  to  allow  the  clouds  to  dwell  for  some 
time  (1-2  seconds)  after  their  initial  formation.  At  the  time  of 
detonation,  the  clouds  had  settled  into  non-hemi spherical  shapes. 
The  dwell  time  that  was  allowed  for  settling  and  evaporation 
resulted  in  explosions  with  relatively  high  effective  blast  yield 
(74  percent  of  the  heat  of  combustion).  However,  it  was  recognized 
that  long  dwell  times  could  not  be  tolerated  in  practice  due  to 
problems  with  cloud  displacement  by  wind.  In  more  recent  tests,  the 
dwell  time  was  reduced.  While  the  clouds  formed  were  not  perfectly 
nemi spherical ,  pressure  waveforms  in  these  tests  were  similar  to 
nuclear  airblast  curve  fits  with  the  exception  that  many  of  the 
waveforms  contained  second  shocks  of  unidentified  origin. 

Early  tests  with  large-scale  jets  from  single  nozzles  achieved 
projection  of  propylene  oxide  to  a  distance  of  about  55  m  (180  ft). 
This  was  short  of  the  71  m  (233  ft)  jet  reach  necessary  to  develop  a 
fuel-air  cloud  at  the  1  KT  nuclear  level.  It  was  felt  that  reach 
was  being  limited  by  the  fuel  pressurization  system  (dispenser) 
which  was  capable  of  pressures  up  to  about  0.5  MPa  (75  psi).  In 
addition,  the  piston  used  in  this  early  dispenser  was  relatively 

small.  As  a  result,  when  large  nozzles  were  used,  the  dispenser 
start-up  time  was  comparable  to  the  total  dispensing  time. 

In  an  attempt  to  achieve  greater  reach,  a  dispenser  having  a 
much  larger  piston  was  designed.  The  piston  in  this  second- 

generation  dispenser  was  driven  by  a  gas  generator  which  provided 

dispensing  pressures  up  to  about  1.6  MPa  (230  psi).  The  start-up 
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time  of  the  new  dispenser  was  sufficiently  fast;  however,  the  anti¬ 
cipated  increases  in  jet  reach  nave  not  as  yet  been  realized.  Based 
upon  experiments  with  smaller  jets,  it  is  currently  believed  that 
while  somewhat  higher  dispensing  pressures  are  still  needed,  the 
primary  obstacle  to  increased  jet  reach  is  the  fuel  capabity  of  the 
di spenser. 

The  present  report  summarizes  progress  made  during  1980  in 
three  areas;  namely,  the  design  of  a  high-capacity,  high-pressure 
impulsive  fuel  dispensing  system,  an  assessment  of  the  possible 
causes  of  secondary  waves  that  have  been  observed  in  fuel -air  explo¬ 
sions,  and  improvements  in  i nstrumentation  which  were  made  to 
enhance  the  quality  of  the  pressure  data.  These  areas  are  discussed 
respectively  in  Sections  II  through  IV.  Concluding  remarks  and 
recommendations  for  further  work  are  given  in  Section  V. 


II.  LARGE  CAPACITY  IMPULSIVE  FUEL  DISPENSER 


A  full-scale  FAE/nuclear  airblast  simulator  at  1  KT  equivalent 
yield  would  impulsively  form  and  then  detonate  hemispherical 
fuel-air  clouds  on  the  order  of  71  m  (233  ft)  radius.  To 
demonstrate  that  individual  fuel  jets  can  be  projected  to  distances 
of  this  magnitude,  a  dispenser  with  considerably  greater  capacity 
than  that  of  the  existing  unit  (Reference  2)  is  required.  A 
dispenser  having  0.42  m3  (110  gal)  capacity  and  14  MPa  (2000  psi) 
operating  pressure  capability  was  designed  and  engineered  during  the 
presently  reported  contract  period.  Details  of  the  design 
considerations  are  discussed  in  this  section. 


1.  BASIC  DESIGN  REQUIREMENTS 

Previous  S  work  with  hydraulic  jets  discharged  from  single 
nozzles  (References  1  and  2)  has  indicated  the  general  influence  on 
jet  reach  of  nozzle  diameter,  jet  discharge  velocity,  and  the 
quantity  of  fuel  dispensed.  Figure  1  illustrates  that,  at  fixed 
nozzle  diameter,  the  maximum  jet  travel  (reach)  first  increases  with 
nozzle  pressure  drop  (or  jet  velocity)  and  then  levels  off  or 
decreases  somewhat.  Above  the  knee  of  the  curve,  which  occurs 
roughly  at  0.6  MPa  (90  psi)  in  Figure  1,  a  substantial  increase  in 
the  degree  of  jet  atomization  and  transverse  jet  spreading  is 
observed.  The  pressure  at  the  knee  is  considered  the  minimum  level 
that  is  useful  for  generating  well  atomized  jets. 

Figure  2  shows  that,  as  would  be  expected,  jet  reach  increases 
with  nozzle  diameter.  Typical  jet  reach  is  in  the  range  of 
1000-2000  nozzle  diameters,  with  1500  being  a  useful  rule-of-thumb 
value  for  hydrocarbon  fuels. 

In  Figure  3,  the  effect  of  the  quantity  of  fuel  dispensed  on 
jet  reach  is  illustrated.  Ihe  jet  reacn  drops  off  rapidly  when  the 
quantity  dispensed  is  below  a  certain  level.  This  is  believed  to  be 
due  to  breakup  of  the  jet  tail  when  the  jet  is  short  (low  dispensed 
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fuel  quantity).  The  disturbance  responsible  for  tail  breakup 
rapidly  propagates  toward  the  jet  nose  (or  head),  causing  premature 
jet  rupture  and  diminished  reach.  A  semi-empirical  relation  for  the 
smallest  volume  of  liquid  that  can  be  dispensed  without  reach 
reduction  was  derived  in  Reference  2.  The  result  is 


7  \  1/2 


Vp  =  0.428  uD 
min 


where  u^  is  jet  discharge  velocity,  is  nozzle  diameter,  and 
,  Op  are  the  fluid  density  and  surface  tension  respectively. 

Although  these  three  basic  features  of  impulsive  jet  breakup 
have  been  investigated  in  survey  experiments,  a  series  of  tests  in 
which  the  controlling  parameters  are  systematically  varied  has  not 
as  yet  been  conducted.  Because  of  this,  the  design  requirements  for 
the  full-scale  single  jet  dispenser  were  necessarily  based  on  some 
rather  crude  estimates. 

A  limited  number  of  experiments  had  earlier  been  carried  out 
with  small  diameter,  steady-state  (continuous)  jets  (Ref.  2).  In 
those  experiments,  it  appeared  that  the  atomization  of  0.36  cm  (0.14 
ini  diameter  propylene  oxide  jets  diminished  substantially  below  jet 
velocities  of  about  25  m/s  (82  ft/s).  For  present  purposes,  it  was 
assumed  that  tnis  velocity  Is  the  minimum  that  would  result  in  full 
reach.  More  recent  tests  with  heptane,  dispensed  through  a  2.5  cm 
(1  in)  diameter  nozzle,  indicated  for  this  size  nozzle  that  the 
minimum  velocity  for  full  reach  is  about  40  m/s  (130  ft/s)  (see 
Figure  1).  Based  on  these  two  data  points,  a  very  rough  extrapola¬ 
tion  to  larget  jet  sizes  was  made.  In  SI  units  this  extrapolation  is 


700  Dn  +  22.5. 


If  jet  reach  is  taken  to  be  1500  nozzle  diameters,  then  a 
nozzle  of  0.06  m  (2.36  in)  diameter  would  be  required  for  a  jet  to 


reach  90  m  (295  ft).  From  Equation  (2),  the  jet  velocity  would  then 
have  to  be  at  least  uQ  =  64.5  m/s  (211  ft/s)  for  the  jet  to  attain 
the  90  m  (295  ft)  reach  and  to  be  at  least  minimally  atomized.  This 
velocity  corresponds  to  a  nozzle  pressure  drop  of  1.42  MPa  (206  psi) 
for  heptane.  Finally,  from  Equation  (1),  at  least  0.264  m  (70 
gal)  of  heptane  (o^  =  0.0209  N/m  (0.00143  lb/ft))  would  have  to  be 
dispensed  through  the  nozzle. 

Due  to  the  considerable  uncertainty  of  these  coarse  estimates, 
it  was  felt  that  the  large  dispenser  should  be  substantially 
over-designed.  For  this  reason  it  was  decided  that  the  dispenser 
would  be  designed  to  operate  at  14  MPa  (2000  psi),  have  a  maximum 

3 

capacity  of  0.42  m  (110  gal),  and  be  capable  of  accommodating 
nozzles  as  large  as  0.1  m  (4  in)  diameter. 

2.  DISPENSER  CONFIGURATION 

Several  possible  dispenser  conf igurations  were  described  in 
Reference  2.  Among  the  various  options,  those  in  which  the  fluid  is 
discharged  by  a  high  pressure  gas  that  is  separated  from  the  fluid 
by  a  piston  were  judged  most  practical.  Working  experience  was 
previously  gained  with  two  dispenser  configurations  of  this  type. 
In  one  (Figure  4),  the  driver  gas  flowed  from  a  source  container 
through  an  orifice  constriction  and  into  a  free  volume  (an  initially 
void  region)  behind  the  piston.  In  concept  this  is  a  good  design 
because  the  pressure  drop-off  after  startup  in  the  free  volume 
(which  increases  during  dispensing)  can  be  minimized  with  a  rela¬ 
tively  small  driver  gas  source  volume.  However,  the  orifice  used 
witn  the  device  in  Figure  4  was  rather  small  so  that  the  initial 
piston  acceleration  was  quite  sluggish.  In  addition,  very  high 
source  pressures  were  needed  to  obtain  the  desired  free  volume 
pressure  levels. 

The  second  dispenser  configuration  tested  employed  a  gas 
generator  to  elevate  the  free  volume  pressure.  This  design  (Figure 
5)  eliminated  the  need  for  a  compressed  gas  source.  An  arrangement 
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Figure  4. 


Schematic  diagram  showing  the  U-tube  used  in  the  experimental 
investigation  for  disseminating  fuel  to  form  hemispherical 
fuel-air  explosive  clouds.  One  leg  of  the  U-tube  was  pres¬ 
surized  in  order  to  force  the  fuel  through  a  nozzle  head 
attached  to  the  other  leg  of  the  U-tube. 
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of  this  type  could  be  quite  successful  if  the  propellant  were 

tailored  to  raise  the  free  volume  pressure  rapidly  at  first,  and  to 
then  burn  more  slowly  so  as  to  sustain  the  pressure  while  dispensing 
progresses.  In  the  experimental  unit  that  was  built,  the  propellant 
used  was  smokeless  powder  (Hercules  Red  Dot).  The  products  of 
combustion  of  this  propellant  were  restrained  in  a  breech  by  a 

diaphragm  until  combustion  was  nearly  completed.  The  products  were 

then  exhausted  impulsively  into  the  free  volume.  This  fast  blowdown 

resulted  in  very  rapid  initial  free  volume  pressure  rise  and  in  good 
piston  acceleration.  However,  during  dispensing  the  free  volume 
pressure  decreased  substantially  due  to  gas  expansion  and  heat 
transfer. 

A  gas  generator  designed  to  meet  the  needs  of  a  full-scale 
simulator  installation  would  appear  at  present  to  be  the  most 
promising  piston  driver  option.  However,  for  experimental 
single-nozzle  work,  it  is  necessary  to  test  a  wide  range  of  nozzle 
sizes,  driver  pressures,  and  dispensed  fuel  quantities.  Since  a 
different  tailored  gas  generator  would  be  required  for  each  such 
variation,  this  approach  would  be  very  costly.  For  that  reason, 
further  use  of  gas  generators  in  the  experimental  work  was  not 
considered  advisable. 

Three  configurations  were  considered  for  the  large  capacity 
dispenser.  These  are  illustrated  schematically  on  Figure  6.  The 
first  (Figure  6a)  would  use  ambient  temperature  air  or  nitrogen  for 
the  driver  gas.  The  high-pressure  gas  source  container  would  be 
directly  linked  to  the  dispenser  cylinder  as  shown.  The  piston 
would  be  initially  restrained  from  motion  by  a  cable  that  would  then 
be  severed  using  an  explosive  actuator  to  start  fluid  dispensing. 
This  configuration  has  the  advantage  that  the  initial  piston 
acceleration  is  maximized  for  any  given  supply  pressure,  since  the 
piston  is  exposed  to  the  full  supply  pressure  in  advance  of  the 
event  that  starts  the  test.  On  the  other  hand,  explosive  actuators 
of  the  size  needed  for  this  application  would  represent  a  nontrivial 


Figure  6.  Large  capacity  dispenser  configurations. 


added  expense.  Additionally,  provisions  would  have  to  be  made  to 
access  tne  area  behind  the  piston  so  as  to  enable  replacement  of  the 
explosive  actuators. 

The  second  dispenser  configuration  that  was  considered.  Figure 
6b,  would  also  make  use  of  a  large-volume  air  or  nitrogen  gas 
supply.  In  this  case,  however,  the  space  (if  any)  behind  the  piston 
would  be  initially  at  ambient  pressure.  A  diaphragm  would  be 
located  in  the  connecting  plumbing  between  the  supply  vessel  and  the 
dispenser  cylinder.  To  start  the  dispensing  process,  this  diaphragm 
would  be  ruptured.  Fast  rupture  would  be  accomplished  through  the 
use  of  a  small  explosive  device,  or  by  means  of  any  of  the 
techniques  used  commonly  in  shock  tubes,  such  as  a  double-diaphragm 
arrangement.  Fast  diaphragm  opening  is  essential  for  adequate 
piston  acceleration  performance;  the  opening  of  an  electrically  or 
air  actuated  valve  is  much  too  slow  for  this  application.  The 
configuration  of  Figure  6b  would  be  simpler  to  build  and  to  operate 
than  that  of  Figure  6a.  However,  the  initial  piston  acceleration 
would  be  somewhat  slower  than  that  in  the  design  of  Figure  6a,  for 
equivalent  driver  pressures.  This  is  due  to  the  time  required  for 
pressure  to  develop  in  the  space  behind  the  piston. 

The  designs  in  Figures  6a  and  6b  both  require  a  considerable 
quantity  of  high  pressure  air  or  nitrogen.  High  pressure  air  could 
be  loaded  with  a  compressor.  The  use  of  such  compressors  requires 
extreme  care,  however,  due  to  the  possibility  for  explosion.  The 
supply  gas  vessel  would  also  have  to  be  kept  extremely  clean  for  the 
same  reason.  In  addition,  a  fairly  large  compressor  would  be  needed 
to  avoid  extended  pump-up  time.  It  is  felt  far  preferable  to  use 
nitrogen  gas  rather  than  compressed  air.  Nitrogen  can  be  supplied 
at  the  time  of  use  from  large  trai lor-mounted  transport  containers. 
In  this  case,  loading  of  the  gas  source  vessel  would  take  only  a  few 
minutes  and  there  would  be  no  hazard  of  explosion.  On  the  other 
hand,  in  the  quantities  required,  the  cost  of  the  nitrogen  and 
delivery  service  would  be  several  hundred  dollars  per  test. 


The  third  dispenser  design  considered  is  that  of  Figure  6c. 
Here  the  arrangement  is  identical  to  that  ot  Figure  6b  except  that 
combustion  products  from  a  gas-phase  fuel-air  explosion  are  used  as 
the  driver  gas.  The  gas  supply  bottle  would  be  charged  with  a 
fuel-air  mixture,  at  relatively  low  initial  pressure,  which  would  be 
detonated  to  start  the  test.  A  diaphragm  would  be  used  (as  in  the 
design  of  Figure  6b);  however,  in  this  case,  only  a  plain,  single 
diaphragm  would  be  needed,  and  the  diaphragm  would  have  only  to 
withstand  a  much  lower  initial  pressure.  Alternately,  the  fuel-air 
mixture  could  be  deflagrated  (ignited  by  a  spark  plug)  and  a  thicker 
diaphragm  employed  which  would  rupture  when  the  combustion  product 
pressure  was  near  its  maximum  level.  The  combustion  driver  scheme 
would  by  far  have  the  lowest  expendible  materials  cost  per  test  of 
the  three  options  in  Figure  6.  It  would  also  have  the  greatest 
safety  if  fuel  injection  were  to  be  done  remotely  just  before  a 
test,  since  there  would  be  no  extended  period  of  high  pressure  gas 
storage.  However,  it  is  not  known  how  accurately  the  fuel  and  air 
can  be  mixed  (if  they  are  loaded  separately);  poor  mixing  would 
compromise  test  repeatibi 1 ity.  Also,  the  supply  gas  pressure  vessel 
would  have  to  be  capable  of  sustaining  the  pressures  developed 
behind  a  reflected  detonation  at  the  vessel  walls. 

To  promote  fuel-air  mixing,  a  number  of  small  fuel  containers 
could  be  located  at  various  points  on  the  pressure  vessel  and 
discharged  into  it  simultaneously.  The  highest  possible  discharge 
velocity  would  be  used  to  encourage  turbulent  mixing  within  the 
vessel.  This  suggests  that  the  fuel  chosen  should  have  a  high  vapor 
pressure  (since  it  would  be  supplied  from  standard  gas  cylinders). 
A  promising  candidate  fuel,  for  example,  would  be  ethylene,  which  is 
relatively  inexpensive  and  has  a  vapor  pressure  or  approximately  8.3 
MPa  (1200  psi). 

The  constant-volume  combustion  characteristics  of  initially 
stoichiometric  ethylene-air  mixtures  were  computed  using  the 
URAKL  code.  The  code  computes  the  product  state  which  is  in 


thermodynamic 

and 

chemical  equilibrium. 

The  result, 
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A  large  pressure  vessel,  located  at  the  S  -operated  Green 

Farm  Test  Site,  was  available  for  use  with  the  high  capacity 

dispenser.  This  is  a  spherical  explosion  vessel  that  had  previously 

been  used  for  confined  detonation  experiments,  and  which  is  capable 

of  withstanding  very  high  pressures.  The  vessel  has  an  internal 
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volume  of  approximately  3.21  m  (113  ft  ).  The  quantity  of 
ethylene  that  it  would  be  necessary  to  inject  into  this  vessel  (with 
the  appropriate  amount  of  air)  was  calculated  from  the  perfect  gas 
law  and  is  plotted  for  reference  on  Figure  8. 

Two  other  simple  calculations  relating  to  the  combustion 
driver  design  are  of  interest.  If  the  fuel-air  mixture  detonates 
within  the  pressure  vessel,  the  pressures  sustained  on  the  vessel 

walls  will  be  considerably  in  excess  of  the  equilibrium,  constant- 
volume  level.  For  stoichiometric  ethylene-air  mixtures,  the 
detonation  pressure  is  about  17.1  times  the  initial  mixture 
pressure.  If  this  detonation  reflects  normally  from  a  surface,  a 
simplified  calculation  shows  that  the  reflected  pressure  could  be 
expected  to  be  about  33.2  times  the  initial  mixture  pressure.  To 

achieve  a  final  equilibrium  pressure  of  13.8  MPa  (2000  psi) 
throughout  the  vessel  requires  an  absolute  initial  ethylene-air 

pressure  of  =  1.44  MPa  (209  psi),  so  that  in  this  extreme  case 
the  reflected  detonation  pressure  would  reach  47.8  MPa  (6930  psi). 
The  spherical  explosion  vessel  mentioned  above  had  been  hydrotested 
when  new  to  52  MPa  (7500  psi)  and  so  it  was  thought  that  it  would  be 
suitable  for  this  service. 

Finally,  an  escimate  of  the  explosion  vessel  wall  temperature, 
after  thermally  equilibrating  with  the  combustion  products,  was  also 
made.  The  explosion  vessel  is  made  of  steel  and  has  a  mass  of 
approximately  1.54  x  10^  KG  (17  tons).  For  the  extreme  case  of 


Equilibrium  product  pressure,  P~,  versus  initial  reactant  mixture  pressure,  P 
constant  volume  combustion  of  ethylene  in  air.  Initial  temperature  is  300°K. 


ethylene  combustion  products  at  13.8  MPa  (2000  psi),  a 
straightforward  First  Law  analysis  gave  an  expected  wall  temperature 
rise  of  approximately  21°C  (38°F)  above  the  initial  ambient  level. 
This  was  regarded  as  negligible. 

All  three  of  the  dispenser  configurations  that  have  been 
described  are  considered  viable.  At  the  time  of  this  writing,  the 
arrangement  of  Figure  6b  is  favored  in  spite  of  its  higher  per-test 
cost.  This  is  because  that  design  appears  to  offer  the  greatest 
experimental  control  with  minimal  operating  complexity.  However,  it 
was  felt  that  the  three  configurations  were  similar  enough  that  the 
design  of  the  dispenser  should  be  such  as  to  accommodate  all  three. 
At  small  additional  cost,  this  would  permit  the  flexibility  to 
change  driver  options  should  such  a  change  become  desirable. 

3.  DISPENSER  SIZING  AND  PERFORMANCE  CALCULATIONS 

The  fundamental  design  variable  in  sizing  the  dispenser  is  the 
diameter  of  the  piston.  For  fixed  dispensed-1 ?quid  capacity  and 
maximum  operating  pressure,  the  cost  of  the  dispenser  assembly 
increases  with  piston  diameter.  As  this  diameter  is  increased,  the 
cylinder  wall  thickness,  the  size  of  the  cylinder  end  caps  and  the 
sizes  of  all  seals  increase.  In  addition,  larger  fabricating 
equipment  is  required,  at  generally  greater  cost.  However,  a  large 
piston  diameter  improves  dispenser  performance  and  minimizes  the 
overall  cylinder  length  (which  simplifies  installation).  Because  of 
these  tradeoffs,  it  was  decided  to  select  the  smallest  piston 
diameter  consistent  with  performance  requirements  and  with  readily 
available  cylinder  lengths. 

The  total  internal  volume  of  the  dispenser  cylinder  is  simply 
the  sum  of  the  volumes  occupied  by  the  piston  and  by  the  maximum 
amount  of  fluid  to  be  dispensed.  In  Figure  9  the  associated  overall 
cylinder  length  is  plotted  as  a  function  of  inside  (or  piston) 
diameter,  assuming  that  the  piston  length  is  twice  its  diameter  (for 
stability).  As  can  be  seen,  rather  unwieldly  cylinder  lengths  would 
be  needed  for  piston  diameters  much  smaller  than  about  30  cm  (12 


Dispenser  performance  was  analyzed  to  first  order  with  a 
simplified  model.  The  sketch  in  Figure  10  defines  the  various 
pressures,  temperatures,  masses,  and  volumes  that  are  involved. 

Since  the  maximum  value  of  Vp  is  much  less  than  Vc,  changes  in 
the  pressure-vessel  pressure  and  temperature  were  ignored. 
Frictional  effects  in  the  connecting  pipe  were  also  neglected  (these 
are  discussed  later  in  this  Section).  The  smallest  area  in  the 
connecting  plumbing,  here  referred  to  as  the  orifice,  is  A  .  The 
fluid  pressure  was  assumed  to  be  uniform  throughout  the  liquid 
reservoir  (the  velocity  throughout  was  taken  to  equal  the  piston 
velocity).  Characteristic  time  intervals  corresponding  to 

significant  changes  in  the  value  of  were  assumed  to  be  long 

compared  to  characteristic  sound  propagation  times  through  the 

liquid.  The  liquid  flow  through  the  nozzle  was  accordingly  treated 
as  quasi-steady,  and  the  Bernoulli  equation  was  used. 

The  equation  of  motion  for  the  piston  under  these  conditions  is 
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where  up  is  the  piston  velocity,  P^  is  ambient  atmospheric 
pressure,  and  is  the  fluid  density.  This  equation  was  derived  in 
Reference  2.  The  rate  of  fluid  efflux  from  the  dispenser  is 
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where  the  discharge  (jet)  velocity  is  (from  mass  conservation) 
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It  is  also  noted  that 
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where 


h  +  VF  "  VT  (7] 

expresses  the  constancy  of  the  total  cylinder  volume  ahead  of  and 
behind  the  piston  at  all  times.  Finally,  the  fluid  pressure  P?  is 
related  to  the  piston  velocity  by 
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from  the  Bernoulli  equation. 

In  order  to  evaluate  Equations  (2)  -  (8),  the  time  varying 
free  volume  pressure,  Pp(t),  must  be  known.  This  was  calculated 
from  thermodynamic  principles,  using  the  expanding  control  volume 
wnich  is  indicated  by  the  dashed  line  on  Figure  10,  and  which 
coincides  with  the  dispenser  free  volume.  For  this  control  volume 
the  energy  equation  is 


dtp  •  O 

dT  =  mi(ei  +  ui/2)  +  miPivi  “  PpupAp 


where  Ep  is  the  total  instantaneous  energy  content  in  the  free 
volume,  and  m.,  e . ,  u^,  P^  and  vi  refer  respectively  to 

the  mass  flow  rate,  specific  internal  energy,  gas  velocity, 

pressure,  and  specific  volume  at  the  free  volume  inlet.  Writing 

Ep  =  rripep  and  noting  that  from  conservation  of  mass  diry/dt  = 

m.,  then  for  a  thermally  and  calorically  perfect  gas,  Equation  (9) 
can  be  reduced  to 


dTr  m.  uAL; 
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where  r  is  the  specific  heats  ratio  of  the  gas.  In  simplifying 
Equation  (9),  quasi-steady  gas  flow  from  the  gas  supply  pressure 
vessel  to  the  free  volume  inlet  was  assumed,  so  that  the  steady 
state  energy  equation  between  those  points,  hc  =  h..  +  u.r/2, 
could  be  used  (h  is  enthalpy).  The  relation  between  Pp  and  Tp 
is  then  simply 

mFRTF 
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where  R  is  the  specific  gas  constant.  Also,  the  equation  for  free 
volume  expansion  is  » 
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To  complete  this  set  of  equations,  the  rate  of  mass  flow  into 
the  free  volume,  m.,  is  needed.  From  mass  conservation,  m.  = 
mo,  where  mo  is  the  flow  rate  through  the  orifice.  If  the 
orifice  is  not  exactly  at  the  free  volume  inlet,  flow  through  it 
will  be  choked  when  Pp  <  Pgg,  where  the  back  pressure  for 
incipient  choking  in  this  case  is 
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and  the  Mach  number  at  the  inlet  corresponding  to  choked  orifice 
flow  is  M^c,  determined  from 


Therefore,  when  Pp 


-  dC’ 
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hand,  if 
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rBC 

Y+l 


choked.  In  this  latter  case,  the  Mach  number  at  the  free  volume 
inlet  is 


! - 

/  pr 

y-l 

M<  -J 

(A) 

k 

Y  -1 

(15) 


and  the  inlet  mass  flow  rate  is 


mi  =pcAiMi  Iwr  ( 1  +  M? ) 


Y+l 

2iY^n 


(16) 


This  group  of  equations  was  numerically  integrated  for  a 

variety  of  possible  conditions.  In  all  cases,  the  initial  piston 

velocity,  up(0)  was  zero,  Pp(0)  =  (0)  =  Pft  =  1  atm,  and  Tp 

=  20°C  (68°F).  The  gas  initially  in  the  free  volume  was  assumed  for 

convenience  to  be  the  same  as  the  driver  gas.  The  initial  mass  in 

the  free  volume  was  calculated  from  Equation  (11).  Required  input 

to  the  calculation  was  the  initial  liquid  and  free  volumes,  (0) 

and  VF(0)  (or  V^);  the  gas  supply  pressure,  temperature,  and 

thermal  properties,  P  ,  T  ,  y,  and  R;  the  fluid  density,  p  •  the 

piston,  nozzle,  orifice  and  free  volume  inlet  cross-sectional  areas, 

Ap,  A^,  Aq,  and  A^ ;  and  the  piston  mass.  The  last  of  these 

was  estimated  from  m  =  p  V  in  which  p  represents  a  bulk 

P  P  P  P 

mass  density  of  the  piston.  The  piston  volume  was  calculated  as 
before  by  assuming  that  the  piston  length-to-diameter  ratio  is  two. 
A  suitable  value  for  p^,  based  on  a  previously  built  dispenser 
piston,  is  pp  =  1150  kg/m3  (71.9  lbm/ft^). 

In  the  calculations,  it  was  assumed  that  the  steel  explosion 
vessel  would  oe  used  as  the  gas  supply  container.  The  largest 
existing  outlet  port  on  this  vessel  is  8.75  cm  (3.44  in)  in 


diameter.  This  size  was  used  in  most  of  the  calculations  as  the 
orifice  diameter.  The  connecting  pipr  between  the  pressure  vessel 
and  the  dispenser  cylinder  inlet  was  assumed  to  be  3-inch 
douole-extra-strong  steel  pipe,  which  has  an  inside  diameter  of  17.5 
cm  (6.37  in).  This  diameter  was  used  to  calculate  the  cylinder 

inlet  area,  A  . 

Figures  11  through  14  illustrate  the  computed  dispenser 

dynamics  for  the  case  of  a  0.102  m  (4  in)  diameter  nozzle  and  a 

u.35o  m  (14  in.)  diameter  piston,  with  V(0)  =  1/3  VT  (all 

3  1 

calculations  assume  =  0.42  m  (110  gal)).  The  gas  supply  in 

tms  case  was  ethylene/air  combustion  products  at  P  =  16.1  MPa 

(  2335  psi),  Tc  =  2876°K  (5177°R),  and  with  =  1.184  and  W  =  28.4 

(molecular  weight) . 

Gas  flow  from  the  supply  vessel  into  the  dispenser  free  volume 
causes  the  free  volume  pressure  Pp  to  begin  rising  at  once.  The 
piston  and  the  fluid  ahead  of  it  begin  to  accelerate,  and  this 
acceleration  initially  increases  with  time,  while  liquid  is  being 

discharged  from  the  dispenser,  the  dispensing  (or  fluid)  pressure  Pj. 
begins  to  rise.  As  this  takes  place  the  piston  acceleration 
steadily  diminishes.  The  piston  velocity  continues  to  increase, 
however,  as  P^,  approaches  but  does  not  exactly  reach  Pp.  During 
this  period,  gas  flow  into  the  dispenser  more  than  compensates  for 
free  volume  expansion,  so  that  Pp  steadily  increases.  This 
process  continues  until  the  gas  flow  to  the  dispenser  ceases  to  be 
choked.  Then  the  inlet  mass  flow  rate  drops  substantially  and 
equilibrium  operation  is  soon  after  established  with  P^  =  Pp  = 
constant  and  up  =  constant  until  all  fluid  has  been  dispensed. 
The  value  of  P^  at  this  time  is  only  slightly  below  the  supply 

pressure.  Note  that  the  temperature  in  the  free  volume  eventually 

increases  above  that  in  the  supply  vessel,  due  to  flow  work.  The 
increase  is  not,  however,  as  great  as  that  which  would  occur  in  the 
absence  of  free  volume  expansion.  (The  free  volume  temperature  is 
of  significance  to  the  selection  of  piston  seals.) 
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wnile  Figures  11  through  14  illustrate  the  dispenser  dynamics 
clearly,  the  example  itself  is  one  of  very  poor  dispenser 


performance.  The  acceleration  or  startup  time  t 


AC’ 


defined  here 


as  the  time  preceding  constant-pressure  dispensing,  is  in  this  case 
nearly  90  percent  of  the  total  dispensing  time,  t^.  This  poor 
performance  is  the  combined  result  of  a  large-diameter  nozzle,  a 
relatively  small  orifice,  and  a  large  initial  free  volume. 

As  the  initial  free  volume  is  reduced,  startup  performance  of 
the  dispenser  rapidly  improves.  The  time  for  the  free  volume  to  be 
rilled  with  supply  gas  decreases  while  the  total  dispensing  time 
increases.  In  Figure  15,  startup  time  is  plotted  as  a  function  of 
the  quantity  of  liquid  dispensed  (see  Equation  (7)).  All  conditions 
other  than  V.  (0)  are  the  same  as  in  the  previous  example.  When  the 
dispenser  is  loaded  to  near  full  liquid  capacity,  the  startup  time 
is  less  than  6  percent  or  the  dispensing  time,  which  is  considered 
satistactory  performance.  On  the  other  hand,  if  the  load  were  much 
less  than  about  1/3  of  full  capacity,  the  dispenser  would  not  reach 
a  constant  pressure  condition  during  the  entire  dispensing  period. 


Although  the  dispenser  was  to  be  designed  for  P ^  =  14  MPa 
(2000  psi),  V  (0)  =  VT  =  0.42  m3  (110  gal),  and  =  0.1  m  (4 
in.),  it  was  anticipated  that  experiments  would  be  necessary  with 
pressures  as  low  as  3.45  MPa  (500  psi)  and  with  fuel  loads  as  low  as 
0.14  m3  (37  gal)  (1/3  capacity).  Since  dispenser  performance 
would  be  poorest  with  1/3-capacity  loads,  it  was  decided  to  further 
analyze  startup  dynamics  only  for  that  case,  so  that  the  dispenser 
could  then  be  designed  for  satisfactory  operation  under  the  most 
difficult  conditions. 

Table  1  summarizes  these  calculations.  In  all  cases,  the 
fluid  dispensed  was  assumed  to  be  heptane.  The  piston  mass  was 
estimated  as  described  earlier.  Th°  niston,  nozzle,  orifice,  and 
dispenser  inlet  diameters,  D  ,  D^,  Dq,  and  Di  are  listed  as 

are  the  driver  gas  pressure  and  temperature,  P  and  T  ,  and  the 
driver  gas  specific  heats  ratio  and  molecular  weight,  y  and  w. 
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Dispenser  performance  is  characterized  principally  by  the  ratio  of 
acceleration  (startup)  time,  t^,  to  total  dispensing  time,  tQ, 
(which  is  also  tabulated  separately  for  reference).  In  all  of  the 
cases  listed,  with  the  exception  of  Case  7,  the  driver  gas  assumed 
was  the  products  of  ethylene-air  combustion. 

Case  Number  1  is  a  reference  base.  Performance  is  quite  poor, 
but  it  is  not  as  bad  as  in  Case  2,  which  was  the  example  discussed 
earlier  (Figures  11  through  14).  The  difference  between  these  two 
cases  is  the  nozzle  diameter,  which  is  7.6  cm  (3  in.)  in  Case  1  and 
10  cm  (4  in. )  in  Case  2. 

Increasing  the  piston  diameter  from  36  cm  (14  in.)  to  41  cm 
(16  in.)  does  not  improve  dispenser  performance  in  this  range  of 
conditions,  as  can  be  seen  by  comparing  Cases  1  and  3.  A  36  cm  (14 
in.)  piston  therefore  appeared  to  be  the  best  design  choice. 

Case  4  demonstrates  the  relative  performance  of  the  dispenser 
at  lower  driver  pressures.  The  ratio  of  startup  to  dispensing  times 
is  considerably  reduced,  compared  with  Case  1.  Although  tftC 
increases  as  Pc  is  reduced,  the  concomitant  increase  in  tQ  is 
much  larger.  Consequently,  worst-case  dispenser  performance  occurs 
at  the  highest  driver  pressures,  contrary  to  what  might  be 
intuitively  expected. 

On  considering  these  first  four  cases,  it  was  evident  that  in 

order  to  attain  acceptable  high-pressure  dispenser  performance,  with 

less  than  full  fluid  capacity,  it  would  be  necessary  to  enlarge  the 

orifice  diameter,  0  .  As  mentioned  earlier,  the  available 

o 

pressure  supply  vessel  has  an  existing  8.74  cm  (3.44  in.)  diameter 
port.  Increasing  the  size  of  this  port  would  be  possible,  but  as 
the  modification  would  represent  additional  expense  and  delay,  it 
had  been  hoped  to  be  able  to  make  use  of  the  existing  opening. 

In  Table  1,  Cases  5  through  7,  the  dispenser  performance  with 
an  enlarged  orifice  is  indicated.  Because  of  the  location  of  the 
explosion  vessel  port,  the  maximum  practical  enlargement  would 


correspond  to  the  inside  diameter  of  standard  8  in.  double-extra¬ 
strong  pipe,  that  is,  17.5  cm  (6.87  in.).  As  can  be  seen  by 

comparing  cases  5  and  1,  such  an  increase  in  Dq  would  greatly 

enhance  dispenser  performance.  The  fractional  startup  time  would  be 
reduced  by  a  factor  of  about  3.5.  This  considerable  improvement 
would  make  dispenser  performance  acceptable  at  1/3  liquid  capacity 

and  the  highest  driver  pressure.  It  is  felt  that  the  pressure 
vessel  modification  to  enlarge  the  port  diameter  would  be  justified. 

Cases  6  and  7  compare  dispenser  performance  at  approximately 
equal  driver  pressure  but  with  different  driver  gases.  In  both 

cases,  the  enlarged  orifice  is  assumed.  Case  6  represents 
performance  with  ethylene-air  combustion  products  at  Pc  =  13  MPa 
(2000  psi).  The  conditions  in  Case  7  are  nearly  identical,  except 
that  the  driver  gas  is  compressed  air  (equilibrated  to  ambient 
temperature).  The  dispenser  performance  using  air  as  the  driver  gas 
is  unfortunately  quite  poor.  The  performance  with  combustion 
products  is  far  superior,  due  to  their  elevated  temperature. 
Although  the  mass  flow  rate  into  the  dispenser  is  lower  with 
combustion  products  than  it  is  with  air,  the  free  volume  pressure 
develops  much  more  rapidly  (see  Equations  (11)  and  (14)).  This 
leads  to  enhanced  piston  acceleration  and  hence  a  shorter 
acceleration  time.  Note  that  the  dispensing  times  are  about  the 
same  in  these  two  cases.  It  is  expected  that  the  performance  of 
nitrogen  as  a  driver  gas  would  be  almost  identical  to  that  of  air. 

4.  FACILITY  LAYOUT 

One  of  the  principal  concerns  in  the  physical  layout  of  the 
large-capacity  dispenser  was  the  ability  to  successfully  photograph 
the  liquid  jets.  In  all  previous  work  the  jets  had  been  discharged 
vertically.  Difficulties  were  encountered  in  photographing  those 
sections  of  the  jets  which  extended  above  the  horizon  in  the 
camera's  field  of  view.  The  jets  could  be  made  visible  on  high¬ 
speed  films  only  when  a  special  film  and  lens  filter  were  conjunc¬ 
tively  used,  and  that  technique  was  successful  only  when  the  sky  was 
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a  deep  blue  (Reference  2).  The  same  jets  were  however  easily 
photographed  without  any  unusual  procedures  when  the  background  was 
dark  green  or  brown  foliage. 

As  a  result  of  these  experiences,  it  was  felt  that  it  would  be 
preferable  to  dispense  the  jets  horizontally.  The  dispenser  would 
be  oriented  so  that  the  fluid  would  be  discharged  in  a  plane 
parallel  to  a  densely  foliated  hillside.  Horizontal  dispensing 
could  be  accomplished  in  two  ways.  Either  the  dispenser  could 
itself  be  mounted  horizontally,  or  the  dispenser  could  be  mounted 
vertically  and  the  jets  directed  horizontally  by  means  of  curved 
nozzles. 

Neither  of  these  options  was  felt  to  be  entirely 
satisfactory.  To  minimize  centrifugal  effects  in  curved  nozzle 
flow,  a  large  radius  of  curvature  would  be  required.  This  would 
result  in  large,  bulky  nozzles.  In  addition,  the  curved  nozzles 
would  have  to  be  rigidly  supported  against  the  substantial  side 
thrust  that  would  be  developed  in  them.  On  the  other  hand,  if  the 
dispenser  were  horizontally  oriented,  the  fluid  would  have  to  be 
initially  confined  by  some  form  of  nozzle  cap.  This  arrangement 
might  prove  satisfactory  if  a  cap  or  diaphragm  were  used  which  would 
definitely  have  no  effect  on  the  jet  flow. 

It  seemed  a  simpler  solution  to  mount  the  dispenser  at  a 
shallow  angle  and  use  straight  nozzles,  as  illustrated  by  the  scale 
drawing  of  Figures  16  and  17.  At  the  location  of  the  spherical 
explosion  vessel,  it  was  possible  to  elevate  the  dispenser  at  a  16 
degree  angle  above  horizontal.  This  was  determined  by  laying  out 
ground  projections  of  the  jet  path  and  perpendicular  camera  line  of 
sight,  using  a  transit.  A  large  helium-filled  balloon  was  then 
raised  to  represent  the  tip  of  a  122  m  (400  ft)  long  jet.  This  is 
about  25  percent  longer  than  the  greatest  actual  jet  reach 
expected.  When  the  balloon  reached  the  top  of  the  background  hill, 
as  viewed  through  the  camera,  its  elevation  angle  with  respect  to 
the  explosion  vessel  was  found  to  be  16  degrees. 
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16.  Dispenser  system  layout  elevation. 


Figure  17.  Dispenser  system  planform. 


When  oriented  at  a  16  degree  elevation  angle,  it  would  be 
possible  to  fill  the  dispenser  to  capacity  without  capping  the 
nozzle  exit.  The  straight  sections  of  the  discharge  nozzles  are  ten 
inside  diameters  in  length.  For  diameters  greater  than  about  3  cm 
(1.2  in.),  the  ends  of  all  nozzles  would  be  above  the  highest  point 
inside  the  dispenser  cylinder. 

On  the  schematic  sketches  of  Figures  16  and  17,  the  dispenser 
cylinder  is  shown  installed  below  grade  adjacent  to  the  existing 
spnerica!  pressure  vessel.  This  arrangement  minimizes  the 
above-grouno  support  structure,  provides  convenient  access  to  the 
nozzle  end  of  the  cylinder,  and  permits  the  reaction  block  to  bear 
against  soil.  The  reaction  block  design  was  not  completed.  However, 
thrust  levels  that  the  block  would  sustain  were  calculated.  In  the 
most  severe  case,  with  Pc  =  13.8  MPa  (2000  psi)  and  =  10.2  cm 
(4  in.),  the  thrust  developed  by  the  dispenser  would  reach  0.225  MN 

4 

1 5 . 07  x  1J  \0). 

An  io.iitionat  concern  relating  to  the  dispenser  layout  was  the 
magnitude  ot  frictional  losses  in  the  connecting  plumbing  between 
the  pres , m  e  n  ,s.- !  and  the  dispenser  cylinder.  Friction  was  not 
considered  ;ti  ’he  performance  analysis  described  earlier.  The 
configuration  sn.wn  in  Figures  16  and  17  would  require  a  minimum  of 
about  6.1  r-i  x2l>  r't!  of  piping  including  three  90  degree  elbows.  The 
largest  pipe  ■  ize  that  can  be  mated  to  the  outlet  port  of  the 
pressure  vtssel  is  8  in.  (double-extra-strong) ,  which  has  an  inside 
diameter  of  17.6  cm  (6.87  in.).  If  necessary,  the  piping  could  be 
enlarged  a  short  distance  beyond  the  outlet  port  to  any  larger 
size.  However,  high  pressure  pipe  is  very  expensive  in  large  sizes 
and  it  was  felt  that  the  smallest  size  that  would  not  significantly 
affect  dispenser  performance  should  be  used. 

A  reasonable  value  for  tne  friction  coefficient  of  common 
steel  pipe  is  f  =  0.0U8.  In  that  case,  each  of  the  three  90  degree 
eloows  contributes  tne  frictional  equivalent  of  about  30  pipe 
diameters  of  length  (Reference  3).  The  total  effective  length  of 


the  connecting  pipe  in  Figures  16  and  17  is  therefore  =  L  + 
90D,  where  L  is  the  actual  length  of  straight  pipe  sections  of 
inside  diameter  D. 

The  pressure  drop  in  the  connecting  pipe  can  be  readily 
calculated  once  the  Mach  number  at  its  entrance  (just  downstream  of 
the  pressure  vessel  port)  is  known.  This  was  estimated  by  ignoring 
the  drop  in  pressure  and  temperature  within  the  pressure  vessel  from 
stagnation  conditions  to  the  pipe  entrance.  The  entrance  Mach 
number  is  then 


where  m  is  the  mass  flow  through  the  pipe  of  area  A,  Pc  and  T 
are  the  driver  gas  stagnation  pressure  and  temperature,  and  Y,  R  are 
the  specific  heats  ratio  and  specific  gas  constant  of  the  driver  gas. 

Once  a  value  for  M^  was  determined,  the  Mach  number  M,,  at 
the  end  of  the  connecting  pipe  (dispenser  cylinder  inlet)  was  found 
from 


Then  the  piping  pressure  loss,  aP  =  P^  -  P ^ ,  was  calculated. 


where  P^  and  P^  are  the  piping  entrance  and  exit  pressures,  and 

P.  =  P„  was  assumed. 

1  c 

Representative  friction  loss  calculations  were  made  for 
dispenser  conditions  corresponding  to  Cases  5  and  7  on  Table  1. 
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These  two  sets  of  conditions  are  similar  except  that  in  Case  5  the 
driver  gas  was  ethylene/air  combustion  products,  while  in  Case  7  the 
driver  gas  was  ambient  temperature  air.  The  mass  flow  rates  used  in 
the  calculations  were  the  values  associated  with  the  constant-pres¬ 
sure  dispensing  period.  The  results  are  plotted  on  Figure  18.  The 
fractional  pressure  drop  due  to  friction  would  be  less  than  one 
percent  for  inside  pipe  diameters  greater  than  D  =  14.2  cm  (5.6 
in.),  if  combustion  products  were  used  as  the  driver  gas.  For 
equivalent  pipe  sizes  the  pressure  losses  are  substantially  greater 
when  air  is  used  as  driver  gas.  This  is  due  to  the  much  larger  mass 
flow  through  the  piping  when  compressed  air  is  used. 

5.  PISTON  BRAKING 

With  a  driver  pressure  of  16.1  MPa  (2335  psi),  the  steady- 
state  velocity  of  a  35.6  cm  (14  in.)  diameter  piston  dispensing 
heptane  through  a  10.1  cm  (4  in.)  diameter  nozzle  is  about  18  m/s 
(59  ft/s).  If  the  piston  mass  were  90  kg  (198  lbm),  its  kinetic 
energy  would  be  14.6  kJ  (1.07  x  10^  ft-lb).  In  order  to  minimize 
the  risk  of  impact  damage  when  the  piston  completes  its  travel,  some 
provision  to  stop  the  piston  in  a  controlled  manner  is  needed. 

The  most  straightforward  braking  device  for  this  application 
is  a  quill-type  shock  absorber,  which  is  illustrated  schematically 
on  Figure  19.  The  quill  is  a  projection  that  attaches  to  the  front 
of  the  piston,  as  shown.  As  the  quill  enters  the  nozzle  inlet,  fluid 
is  forced  by  the  advancing  piston  to  flow  from  the  reservoir 
indicated  through  a  continually  decreasing  minimum  area.  If  the 
quill  length  and  shape  are  appropriately  designed,  the  reservoir 
pressure  will  rise  during  this  process  in  such  a  way  that  the  piston 
is  smoothly  decelerated  to  zero  velocity. 

During  the  period  of  steady-state  dispensing,  the  piston 
velocity  is  constant  and  the  pressures  acting  on  it  are  uniform  and 
equal  on  both  ends.  As  the  quill  enters  the  nozzle  inlet,  the 
pressure  on  the  back  side  of  the  piston  remains  uniform  and  at  the 
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Figure  19.  Schematic  of  quill  shock  absorber  piston  braking  device. 


level  Pp  (see  Figure  10).  During  this  period  the  pressure  in  the 
reservoir,  Pp,  rises,  while  the  fluid  pressure  drops  off  along  the 
quill  surface  to  a  value  below  Pp  at  its  front  end.  The  reduced 
pressure  on  the  front  surface  of  the  quill  opposes  piston 
deceleration,  so  that  it  is  desirable  to  minimize  the  area  of  that 
surface.  In  the  simplified  analysis  that  follows,  the  pressure  on 
the  entire  surface  of  the  quil1  was  assumed,  for  purposes  of  writing 
the  piston  equation  of  motion,  to  be  zero.  This  assumption  was 
expected  to  result  in  calculated  reservoir  pressures  somewhat  above 
those  actually  needed. 

Accordingly,  the  piston  motion  during  the  deceleration  period 
is  described  by 


mp  dt 


PCA  - 
F  P 


PRAA 


(20 


and 


_ 

dt 


(21 


in  which  A^  is  the  annular  area  between  the  base  of  the  quill  and 
the  inside  wall  of  the  cylinder  body.  The  time  t  =  0  is  defined  as 
the  moment  of  first  entry  by  the  quill  into  the  nozzle  inlet,  which 
is  the  condition  shown  in  Figure  19.  At  this  moment  the  piston 


position,  x 


is  taken  to  be  zero,  and  the  piston  velocity  is  u 


=  u„ 


The  mode  of  stopping  selected  was  that  of  constant 
deceleration.  This  would  require  Pp(t)  =  constant  in  Equation 
20.  In  reality,  Pp  could  be  constant  only  after  a  short  initial 
interval  during  which  it  increases  from  Pp,  but  this  initial 
interval  was  ignored,  equation  (20)  then  integrates  to 


u 


P 


Bt 


(22) 
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where 


praa  pfap 


The  piston  will  have  stopped  (Up  =  0)  at  the  time 


tE  =  B 


From  Equations  (21)  and  (22)  the  piston  path  is  simply 


1  2 

x  =  u  t  -  4  Bt4 
P  PQ  2 


so  that  at  t  =  tr,  the  piston  will  have  moved  the  distance 


Now,  equating  x  in  Equation  (24)  to  the  length  of  the  quill, 

PE 

Lq,  and  using  Equation  (23),  leads  to  the  required  reservoir 
pressure, 

m  uP  n  . 

P  Pn  PFAn 

p _ °  +  t-E  (25) 

R  2LqAa  Aa  '  ' 

The  required  surface  geometry  of  the  quill  was  established  by 
considering  the  flow  out  of  the  reservoir.  It  was  assumed  that  the 
pressure  is  zero  at  and  downstream  of  the  minimum  flow  area,  A^, 
which  occurs  between  the  quill  and  nozzle  inlet  surfaces.  If  it  is 
further  assumed  that  the  flow  is  quasi-steady,  then  the 
instantaneous  mass  flow  rate  is 
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However,  if  the  fluid  is  considered  incompressible,  the  mass  flowing 
through  the  minimum  area  must  equal  that  displaced  by  the  moving 
piston;  that  is, 


mc  =  Vpaa  ■  (27 

Equating  Equations  (26)  and  (27)  to  solve  for  A^,  and 
incorporating  Equation  (25),  gives 

am  =  pPUDAA  / - P -  •  (28 

'pA<vp/VfV 

The  piston  velocity  is  the  only  time-dependent  variable  in 
Equation  (28).  It  is  related  to  its  instantaneous  position  by 

du„  du„  .. 

p  _  p  dt 

dx  dt  dx  ’ 

P  P 

or,  in  view  of  Equations  (20)  and  (21)  , 


du„  PCA  -  PdA. 
p  _  F  p  R  A 

dx  mu 

P  P  P 


Upon  integrating,  and  using  Equation  (25),  this  becomes 


V  up0  (1  -  VV; 


Therefore  Equation  (28)  can  be  written 


This  equation  is  an  explicit  relation  between  the  minimum  flow 
area,  A^,  and  the  piston  position,  x  ,  during  the  constant 

deceleration  process.  It  is  clearly  impossible  with  a  single, 
fixed-surface  quill  to  accomplish  uniform  stopping  under  all 
dispenser  operating  conditions.  The  quill  was  therefore  designed  to 
be  replaceable. 

The  calculated  variation  in  with  xp  appears  in  Figure 

20  for  the  following  case:  Lq  =  0.152  m  (6  in.),  A^  =  0.0813 

m2  (126  in2),  =  683  kg/m3  (42.7  lbm/ft3),  mp  =  91  kg 

(200  lbm) ,  up  =  18  m/s  (59.0  ft/s),  Ap  =  0.0994  m2  (154 

in2),  and  Pp  °=  13.8  MPa  (2000  psi).  The  areas  Aft  and  Ap 

correspond  to  a  piston  diameter  of  35.6  cm  (14  in.)  and  a  quill  base 
diameter  of  15.2  cm  (6  in.).  The  liquid  density  is  that  of 

heptane.  Under  these  conditions  the  reservoir  pressure,  calculated 
from  Equation  (25),  would  be  P^  =  18.1  MPa  (2620  psi)  and  the 
braking  time  would  be  tp  =  16.9  ms.  The  variation  in  quill  radius 

along  its  length  was  determined  from  the  data  in  Figure  20  by 

graphical  means. 


Figure  20.  Minimum  flow  area  between  quill  and  nozzle  inlet  as  a  function  of  piston  travel  during  decel¬ 
eration  period  (see  text  for  conditions). 


III.  SECONDARY  WAVE  SYSTEMS  IN  FUEL-AIR  EXPLOSIONS 
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One  of  the  issues  of  FAE/nuclear  airblast  simulation 
feasibility  that  has  not  been  fully  resolved  is  the  question  of 
pressure  waveform  fidelity.  By  this  is  meant  the  degree  of  temporal 
correlation  between  the  static  pressures  produced  at  equivalent 
radii  in  FAE  and  in  scaled  nuclear  airblasts.  Accurate  waveform 
correspondence  is  considered  requisite  for  meaningful  airblast 
testing  of  certain  targets. 

Four  types  of  waveform  aberrations  are  illustrated 
schematically  on  Figure  21.  Negative  curvature  in  the  pressure 
decay.  Figure  21(a),  leads  to  early  impulse  buildup  as  compared  with 
positive  curvature  (shown  as  a  dashed  line).  The  spiked  excursion 
shown  in  Figure  21(b)  would  have  a  slapping  effect  on  a  target,  but 
it  does  not  otherwise  significantly  alter  the  average  pressure 
profile  or  waveform  impulse.  The  low-amplitude  secondary  wave  in 
Figure  21(c)  occurs  in  the  negative  phase  and  is  not  normally 
considered  deleterious.  The  source  of  this  particular  wave  is  well 
understood  (e.g.,  see  Reference  2). 

The  fourth  type  of  waveform  distortion,  represented  by  Figure 
21(d),  is  presently  of  greatest  concern.  Pressure  excursions  of 
this  type,  which  occur  during  the  positive  phase,  consist  either  of 
compression  wave  packets  or  of  shock  fronts,  followed  by  relatively 
slow  expansion.  These  large-amplitude  secondary  wave  systems 
contribute  non-negl igible  additional  overpressure  impulse  to  the 
blast. 

A  particularly  dramatic  example  of  a  large-amplitude  secondary 
wave  is  given  in  Figure  22.  Secondary  waves  of  this  kind  have  been 
observed  in  many  of  the  FAE  simulator  tests.  Most,  however,  are  not 
of  the  magnitude  in  Figure  22,  and  it  has  frequently  been  difficult 
to  distinguish  what  are  thought  to  be  secondary  waves  from  other 
waveform  disturbances. 
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Types  of  pressure-time  waveform  aberrations  in  fuel-air  e>plosions  (schematics);  (a)  negative 
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Figure  22.  Example  of  large  amplitude  secondary  wave.  Arrow  indicates  conjectured  location  of  wave 
front.  Data  from  FAE  test  CD-14;  22.7  kg  (50  lbm)  propylene  oxide,  R  =  10.1  m  (33.0  ft) 


The  static  pressure  records  in  Figures  23  to  28  are  unsmoothed 
raw  data  taken  from  a  recent  test  of  the  nominal  1/4-tone  FAE  simu¬ 
lator.  In  this  test  it  appears  that  a  compression  wave  packet, 
originating  within  the  exploded  fuel-air  cloud,  coalesced  into  a 
shock  while  overtaking  the  blastwave  front.  An  x-t  plot  of  the 
secondary  wave  front  appears  in  Figure  28  (what  is  believed  to  be 
the  front  is  indicated  by  arrows  on  Figures  23  to  27).  A  much 
greater  spatial  resolution  of  blastwave  pressures  is  needed  in  order 
to  identify  and  trace  the  progress  of  secondary  waves  with  greater 
accuracy.  For  present  purposes,  secondary  wave  fronts  have  been 
associated  with  those  points  on  the  pressure  records  at  which  abrupt 
increases  in  either  the  average  pressure  levels  or  the  average 
values  of  dP/dt  are  observed. 

The  secondary  waves  appear  to  be  somewhat  rotational ly  asym¬ 
metrical.  Figure  30  is  a  pressure  record  from  the  same  test  and 
approximately  the  same  radial  distance  as  that  of  Figure  25,  but  the 
two  pressure  transducers  were  90  degrees  apart.  On  Figure  30,  the 
front  of  the  secondary  wave  seems  relatively  clear,  at  least  as 
compared  with  Figure  25.  However,  when  the  two  pressure  records 
were  superimposed.  Figure  31,  what  had  appeared  to  be  the  respective 
secondary  wave  fronts  did  not  align.  There  does  seem  to  be  rea¬ 
soned  e  rotational  symmetry  in  blastwave  time  of  arrival,  peak 
overpressure,  positive  phase  duration,  and  impulse  development 
(Figure  32 j.  There  also  appears  to  be  some  degree  of  pressure-time 
correspondence  over  that  segment  of  the  secondary  wave  marked  by 
arrows  on  Figure  31. 

In  most  cases,  static  and  stagnation  pressure  records,  taken 
at  the  same  location,  are  consistent  with  one  another,  indicating 
tne  arrival  and  passage  of  secondary  waves  at  nearly  coincident 
times.  An  example  if  Figure  33.  However,  in  some  instances,  such 
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Figure  23.  Secondary  wave  in  experimental  pressure  profile.  Arrow  indicates  approximate  location  of 
wave  front.  Oata  from  FAE  test  CO-13;  22.7  kg  (50  lbm)  propylene  oxide,  R  =  2.91  m  (9.54 
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Figure  30.  Secondary  wave  in  experimental  pressure  profile.  Arrow  indicates  conjectured  location  of 
wave  front.  Data  from  FAE  test  CD-13;  22.7  kg  (50  lbm)  propylene  oxide,  R  =  10.1  m  (33.0 
Transducer  located  90°  from  that  of  Figure  25. 
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Figure  31.  Superposition  of  Figures  25  and  30.  Time  referred  for  both  traces  to  same  instant 
(detonation  initiation). 


as  in  Figure  34,  the  two  waveforms  do  not  correspond,  for  reasons 
that  are  not  presently  understood. 

A  possible  source  of  the  large-ampl itude  secondary  wave 
systems  could  be  the  result  of  temporally  or  spatially  incoherent 
energy  release  by  the  detonating  fuel -air  cloud.  The  energy  of 
combustion  is  apparently  being  coupled  to  the  atmosphere  asym¬ 
metrically,  or  as  a  series  of  pulses,  rather  than  added  uniformly 
and  only  behind  a  point-symmetrical  detonation  front  of  zero  thick¬ 
ness.  (Calculations  of  ideal  fuel-air  explosions  (References  2,  4, 
5)  have  generated  only  low-amplitude,  negative-phase  waves  of  the 
type  represented  schematically  on  Figure  21(c).  All  of  these  calcu¬ 
lations  have  in  common  the  assumptions  of  point  symmetry;  uniform 
initial  fuel -air  mixture;  energy  release  only  at  the  detonation 
front;  kinetically  frozen  combustion  produces;  perfect  gas  behavior; 
and  central,  hard  detonation  initiation.  Large-amplitude  secondary 
waves  have  also  not  been  observed  during  FAE  and  FOE*  experiments  in 
wnich  the  clouds  were  initially  envelope-confined  (References  6 
through  8),  with  one  possible  exception  (Reference  9).  The  use  of 
envelope  confinement  in  these  experiments  assured  accurate  geometric 
symmetry  and  precise  central  initiation. 

Several  characteristics  of  fuel-air  explosions  are  thought  to 
have  a  possible  influence  on  the  formation  of  secondary  waves. 
These  are  cloud  shape,  initiator  location,  fuel  distribution, 
initiator  strength,  detonation  reaction  zone  length,  and  product  gas 
chemical  equilibrium  shifts. 

1.  CLOUD  SHAPE 

If  the  detonating  fuel-air  cloud  is  not  hemispherical,  regions 
adjacent  to  the  cloud  boundaries  will  exist  in  which  components  of 
the  detonation  and  the  transmitted  air  shock  meet  at  oblique  angles, 
due  to  differences  in  their  propagation  velocities  and  directions. 
The  oblique  shock  interactions  will  generate  reflected  waves  which 


*FueT-oxygen  explosion. 


may  be  shocks.  The  reflected  waves  could  strike  the  ground  at  an 
appreciable  distance  behind  the  lead  shock  if  the  interaction  occurs 
substantially  above  ground. 

This  phenomenon,  due  to  shape,  has  been  observed  with 
cylindrical  HE  charges  having  sharp  corners  (References  10  and  11). 
The  mechanism  is  illustrated  schematically  for  that  case  on  Figure 
35,  while  Figure  36  shows  a  series  of  experimental  pressure  records 
taken  on  the  ground  along  the  axis  of  a  detonated  cylindrical  HE 
charge.  A  strong  secondary  wave  can  be  seen  on  the  first  three 
traces  as  it  overtakes  the  leading  blastwave  shock.  According  to 
Reference  10,  such  secondary  waves  were  eliminated  when  the  ends  of 
the  cylindrical  charges  were  shaped  into  hemispheres.  This 
demonstrated  that  the  sharp-cornered  charge  shape  was  responsible 
for  the  secondary  waves.  It  is  possible  that  an  equivalent 
mechanism  could  be  operative  in  fuel -air  cloud  explosions,  if  the 
clouds  formed  are  not  exactly  hemispherical. 

2.  INITIATOR  LOCATION 

Ideally,  an  FAE  cloud  should  be  initiated  at  the  center  of 
symmetry.  If  a  normal  hemisphere  is  initiated  above  g  jund,  either 
a  spherical  detonation  or  the  initiating  blastwave  will  reflect  from 
the  ground.  The  region  of  reflection  will  comprise  an  oblique 
detonation  front  and  an  oblique  shock  front.  This  differs  from  the 
ordinary  reflection  of  a  spherical  blastwave  from  a  surface,  and  may 
be  accompanied  by  additional  shocks  or  expansions.  In  tt. '  case  of 
marginal  detonation  initiation,  partially  reacted  or  unreacted  fuel 
may  be  ignited  by  the  reflected  shock,  producing  a  second  reaction 
front.  In  addition,  above-ground  initiation  would  result  in  a 
spherical  detonation  that  would  overtake  the  cloud  boundary  in  a 
non-coinc ident  fashion.  This  could  produce  a  boundary  interaction, 
leading  to  reflected  waves,  in  a  manner  equivalent  to  that  described 
in  (1)  above.  The  same  type  of  boundary  interaction  would  also 
occur  if  the  initiator  were  located  on  the  ground  but  were  offset 
horizontally.  Horizontal  misalignment  of  the  initiator  is 
equivalent  to  cloud  displacement  due  to  atmospheric  movement. 
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Figure  36.  Experimental  pressure  records  from  detonation  of  cylindrical 
explosive  charge  with  sharp  corners.  Data  taken  on  ground 

along  charge  axis  direction  (from  Reference  11).  * 
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Figure  37  is  a  frame  from  a  high  speed  film  of  a  propylene 
oxide  FAE.  The  25  gm  (0.55  lbm)  C-4  initiator  in  this  test  was 
located  1.2  m  (4  ft)  above  ground  along  the  vertical  axis  of  the 
cloud.  The  reflected  wave  is  very  clearly  defined  in  this  example. 
It  should  be  noted  that  the  strength  of  the  reflected  wave  is 
associated  not  only  with  the  initiator  energy,  which  was  very  small 
in  tnis  case,  but  with  all  of  the  detonated  fuel  in  the  sphere  of 
radius  equal  to  the  initiator  height  (vertical  offset). 

3.  FUEL  DISTRIBUTION 

Non-uniform  fuel  distribution  would  lead  to  refraction  of  the 
detonation  front  within  the  cloud.  This  could  generate  reflected 
compression  or  expansion  waves.  If  the  initiator  strength  were 
marginal,  initiation  might  be  directional,  leading  to  an  irregular 
detonation  front  with  multiple  wave  reflections. 

Fuel  non-uniformity  could  take  several  forms.  The  fuel  could 
be  contained  in  pockets  spaced  irregularly  within  the  cloud.  This 
type  of  non-uniformity  could  account  for  rotational  asymmetry  of 
secondary  waves.  Another  non-uniformity  might  be  a  large  void  volume 
in  the  vicinity  of  the  cloud  center,  which  would  lead  to  multiple 
wave  reflections  (presuming  that  the  cloud  could  be  centrally 

initiated  under  these  conditions).  Also,  the  droplet  size 

distribution  might  vary  within  the  cloud,  so  that  the  detonation 
reaction  zone  would  be  continually  undergoing  transitions  with 

accompanying  wave  reflections. 

4.  INITIATOR  STRENGTH 

Critical  region  (marginal)  initiation  (Reference  1?)  could 
potentially  generate  one  or  more  nested  explosions,  as  well  as  an 
inward-moving  detonation  (retonation) ,  similar  to  what  occurs  in 

tief lagration-to-detonation  transition  (Reference  13).  When 
initiation  is  marginal,  a  galloping  behavior  is  observed  (Reference 
14);  that  is,  the  detonation  front  velocity  oscillates  in  a  pulsed 


Figure  37.  Frame  from  high  speed  film  taken  during  FAE  test  2533; 

22.7  kg  (50  Ibm)  propylene  oxide.  Initiator  height  was 
1.22  m  (4  ft)  above  ground,  located  on  cloud  axis.  The 
distance  between  the  vertical  poles  is  9.2  m  (30  ft). 


ed 


manner.  Presumably,  reflected  waves  are  formed  in  this  process.  In 
addition,  in  the  case  of  marginal  initiation,  the  detonation  shock 
front  and  reaction  zone  tend  to  separate.  If  the  reactions 
following  ignition  are  then  extremely  fast,  a  compression  pulse  will 
form  during  this  unsteady  initiation  process. 

5.  REACTION  ZONE  LENGTH 

If  initiation  is  not  margina,  or  excessive,  a  steady 
detonation  develops  within  the  boundaries  of  the  FAE  cloud.  The 
normal  hydrodynamic  structure  of  a  detonation  comprises  a  shock 
front  followed  by  a  pressure  drop  when  the  smooth,  exothermic 
cna i n-break i ng  reactions  commence.  However,  in  the  case  of  some 
1  iquia-fuel /oxidizer  combinations  (especially  when  liquid  fuels  and 
pure  oxygen  are  involved)  the  ignition  is  not  smooth  but  explosive 
(Reference  15).  In  such  cases,  multiply-nested  blastwaves, 
superimposed  on  a  globally  decaying  pressure  profile,  appear  in  the 
reaction  zone  (Reference  16).  If  the  reaction  zone  is  extremely 
long  due  to  large  droplets,  marginal  stoichiometry,  or  substantial 
out  not  complete  pre-vapor ization  of  the  fuel  (Reference  17),  these 
nested  blastwaves  could  persist  beyond  the  cloud  boundary. 

A  second  effect  of  a  long  reaction  zone  might  be  that  it  would 
produce  a  strengthened  rarefaction  behind  the  detonation.  A 
strengthened  rarefaction  could  be  expected  if  the  mass  in  the 
detonation  reaction  zone  is  comparable  to  the  total  fuel-air  cloud 
mass  ^Reference  2).  A  strong  rarefaction  would  change  the 
characteristics  of  the  ^as  dynamics  which  normally  generate  an 
inconsequential  secondary  wave  in  the  negative  phase.  It  may  be 
that  this  secondary  wave  would  be  amplified  by  the  altered  gas 
dynamics. 

h.  C0M6USTIUN  PRODUCT  EQUILIBRIA 

The  products  of  combustion  in  a  detonation  are  highly 
dissociated  due  to  their  elevated  temperature.  In  the  rarefaction 


that  follows,  the  product  temperature  drops  substantially.  Further 
temperature  reduction  occurs  within  a  fuel-air  cloud  after  the 
detonation  passes  through  the  cloud/air  interface.  Equilibrium 
recombination  of  the  combustion  products  accompanies  these 
temperature  reductions,  with  recovery  of  the  dissociation  energies. 
If  the  recombination  occurs  with  sufficient  rapidity,  a  compression 
wave  packet  could  be  formed  from  this  effectively  impulsive  energy 
addition. 

As  an  example.  Figure  38  shows  the  calculated  time  variation 
of  static  temperature  within  a  detonated  heptane  vapor-air  cloud. 
The  calculation  (Reference  2)  was  carried  out  for  a  uniformly 
stoichiometric,  4.42  m  (14.5  ft)  radius  cloud,  which  corresponds  to 
tne  small  scale  experiments.  The  products  of  combustion  in  this 
calculation  were  assumed  chemically  frozen  so  that  the  temperature 
variations  indicated  are  only  suggestive  of  the  magnitude  of 
change.  The  detonation  front  arrives  at  the  location  in  Figure  38, 
R  =  2.91  m  (9.54  ft),  at  t  =  1.62  ms.  At  that  time,  the  temperature 
rises  to  the  C-J  detonation  level,  T  =  2550  °C  (4620°F).  The 
detonation  products  consist  of  significant  amounts  of  the  chemical 
species  N^,  H^O,  CO,  CO^,  0^,  OH,  H^,  NO,  H,  and  0. 
Within  the  next  8  to  10  ms  the  temperature  at  this  location  drops  by 
nearly  1000°C  (1800°F)  to  about  T  =  1580°C  (2870°F).  At  this 
relatively  low  final  temperature  the  products  will  re-equilibrate  to 
a  condition  in  which  the  major  constituents  are  N^,  H^O,  and 
CC^.  This  represents  effectively  complete  combustion  and  the 
energy  of  dissociation  expended  at  the  higher  temperature  is 
returned  to  the  flow. 

It  was  anticipated  that  the  effects  of  initiator  location, 
initiator  strength,  equilibrium  shifts,  and  detonation  reaction  zone 
length  would  likely  be  negligible  with  large-scale  clouds,  although 
these  same  effects  might  be  contributory  or  dominant  at  smaller 
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Figure  38.  Calculated  temperature  variation  inside  detonated  fuel-air 
cloud  (see  text  for  conditions). 


scales.  However,  it  was  felt  that  any  impact  of  cloud  shape 
asymmetry  or  fuel  distribution  non-uniformity  would  not  diminish 
with  increasing  scale,  and  for  that  reason  it  was  thought  that  their 
effects  should  be  investigated  first. 

The  shape  of  the  fuel-air  cloud  is  important  to  the  nuclear 
airolast  simulator  objectives  for  reasons  other  than  the  possible 
contribution  of  asymmetry  to  secondary  waves.  A  hemispherical  FAE 
cloud  is  theoretically  the  most  efficient  shape  (Reference  13)  and 
the  cloud  dimensions  are  minimized  for  a  given  specified  yield.  In 
addition,  predictable  testing  in  all  directions  and  at  the  highest 
overpressure  levels  is  possible  from  a  hemispherical  FAE.  However, 
truly  hemispherical  clouds  have  not  as  yet  been  formed  at  small 
scale  witn  the  point  source,  multiple-nozzle  impulsive  dispenser. 
It  was  accordingly  felt  that  the  initial  attempts  to  eliminate 
secondary  waves  in  FAE  airblasts  should  be  directed  to  improvements 
in  cloud  shaping.  If  cloud  shape  were  to  prove  inconsequential  to 
the  formation  of  secondary  waves,  the  FAE  simulator  technology  would 
in  any  case  be  advanced  by  such  investigations. 

The  shape  of  FAE  clouds  that  are  formed  by  point-source 
dissemination  is  distorted  by  large  scale  air  vortices  (Reference 
2j.  The  process  by  which  this  distortion  takes  place  is  illustrated 
on  Figure  39.  Liquid  jets  flowing  radially  outward  from  the  source 
disseminator  drag  on  the  atmospheric  air,  causing  it  to  move  away 
trom  the  source.  This  reduces  the  pressure  near  the  center  of  the 
cloud  so  that  air  is  drawn  inward  on  the  underside  of  it.  At  the 
center,  this  inward-moving  air  is  turned  vertically,  forming  into 
torroidal  vortices  around  the  cloud  axis.  Droplets  become  entrained 
in  the  vortex  motion,  deforming  the  cloud.  The  vertical  air 
movements  in  the  vicinity  of  the  dispenser  deflect  the  liquid  jets 
as  they  emerge  from  the  nozzles,  further  distorting  the  cloud 
snape.  An  example  of  a  distorted  cloud  is  given  in  Figure  40. 


Cloud  shape  distortion  is  aggravated  by  increasing  either  the 
dispenser  jet  momentum  or  the  number  of  jets  in  the  hemispherical 


ortions  of  expanding  point-source  spray  clouds  induced  by  aerodynamic  interact 
ipp n  tinrav  and  air.  The  dashed  lines  represent  the  desired  hemispherical  shape 


Frame  from  high  speed  film  showing  fuel -air  cloud  distorted 
by  large-scale  air  vortices;  test  CF-9,  18.1  kg  (39.9  lbm) 
heptane,  1400-nozzle  point  source.  The  distance  between 
the  diamond-shaped  markers  is  9.2  m  (30  ft). 
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array.  When  a  small  number  of  low  momentum  jets  are  used, 
large-scale  air  vortices  do  not  form  because  the  smaller  pressure 
drop  at  the  cloud  center  is  relieved  by  air  counterflow  between 
jets.  In  this  case  it  is  likely  that  a  large  number  of  low-strength 
vortices  develop  between  pairs  of  adjacent  jets.  Unfortunately,  the 
fuel  is  poorly  dispersed  with  this  arrangement,  being  confined  to 
cylindrical  regions  near  the  individual  jet  centerlines.  The  cloud 
takes  on  a  spiked  appearance,  an  example  of  which  appears  in  Figure 
41.  When  detonated,  clouds  of  this  shape  generate  explosions  of 
relatively  low  efficiency.  On  the  other  hand,  if  a  greater  number 
of  liquid  jets  at  higher  velocity  are  used,  the  jets  spread 
laterally  and  overlap.  This  fills  in  the  cloud  more  uniformly,  but 
it  is  when  the  cloud  is  filled  in  that  large-scale  voctices  form  and 
distort  the  cloud  shape  as  in  Figure  40. 

A  number  of  attempts  were  made  to  overcome  vortex-induced 
cloud  distortion  while  maintaining  reasonably  uniform  fuel 
distribution.  For  these  tests,  a  large  number  of  high  velocity  jets 
was  used,  but  several  clusters  of  dispenser  nozzles  were  blocked. 
The  blocked  groups  of  nozzles  were  uniformly  spaced  over  the 
hemispherical  array.  The  purpose  of  blocking  clusters  of  nozzles 
was  to  provide  a  limited  number  of  unobstructed  air  paths  to  the 
cloud  center.  It  was  hoped  that  a  blockage  pattern  could  be  found 
that  would  allow  enough  airflow  through  these  openings  as  to  relieve 
the  reduced  pressure  at  the  cloud  center,  without  compromising  the 
fuel  distribution  uniformity  throughout  the  majority  of  the 
remaining  cloud  volume.  This  would  reduce  the  tendency  for  air  to 
spontaneously  flow  into  the  cloud  on  its  underside,  and  would 
thereby  prevent  the  formation  of  large-scale  torroidal  vortices. 

The  two  types  of  nozzle  blockage  patterns  that  were  tested  are 
illustrated  on  Figure  42.  Blockage  was  accomplished  by  bolting 
machine  screws  through  the  nozzle  openings  in  the  dispenser  head. 
The  dimple  pattern  in  Figure  41(a)  consists  of  seven  equally  spaced 
circular  clusters.  Tests  were  conducted  with  this  pattern  in  which 
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F,  ime  from  high  speed  film  showing  fuel -air  cloud  formed 
from  partially  overlapping  jets;  FAE  test  CD-9,  78.5  percent 
heptane,  21.5  percent  propyl  nitrate  by  volume,  0.0265  m3 
(7  gal.)  total  mixture,  600-nozzle  point  source.  The  dis¬ 
tance  between  the  diamond-shaped  markers  is  9.2  m  (30  ft). 
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Figure  42.  Dispenser  head  nozzle  blockage  patterns:  (a)  dimple-type 
pattern,  (b)  ring  pattern. 


the  total  blockage  amounted  to  2,  10,  15,  20,  and  25  percent  of  the 
total  number  of  nozzles.  The  other  blockage  pattern,  which  is  shown 
in  Figure  42(b),  is  a  continuous  axisymmetric  ring.  The  ring  was 
located  so  as  to  roughly  divide  the  nozzle  head  into  two  equal  areas 
(with  the  same  number  of  nozzles  above  and  below  the  ring).  In  this 
pattern,  tests  were  conducted  with  blockage  amounting  to  5.5,  7.5, 
and  10  percent  of  the  total  number  of  nozzles.  The  test  conditions 
are  summarized  in  Appendix  B. 

Figures  43  and  44  show  the  effect  on  cloud  shape  of  nozzle 
blockage  in  the  pattern  of  Figure  42(a).  These  examples  should  be 
compared  with  Figure  40,  which  shows  a  cloud  formed  under  similar 
conditions  but  with  no  blockage.  With  25  percent  blockage  (Figure 
43),  the  cloud  extends  fully  to  the  ground,  and  its  overall  size  is 
greatly  increased.  The  shape  is  characterized  by  a  number  of  large 
thumblike  projections  and  deep  indentations,  presumably  the  result 
of  local  vortices  that  form  from  air  drawn  in  through  the  blocked 
areas.  When  a  cloud  of  this  shape  was  detonated,  the  pressure 
records  on  Figure  45  were  obtained.  A  secondary  wave  of  moderate 
strength  appears  on  these  records.  Reducing  the  nozzle  blockage  to 
10  percent  seemed  to  improve  the  cloud  shape,  but  when  the  cloud  in 
Figure  44  was  detonated,  the  pressure  records  (Figure  46)  indicated 
the  presence  of  a  strong  secondary  wave. 

Nozzle  blockage  in  the  ring  pattern  (Figure  42(b))  yielded 
very  striking  cloud  shape  changes.  The  changes,  however,  were  not 
in  the  direction  of  improved  hemisphericity.  Figure  47  is  an 
example  of  this.  The  cloud  extends  to  the  ground  but  the  void 
region  is  very  large  compared  with  the  fractional  area  that  was 
blocked  on  the  nozzle  head.  None  of  the  clouds  having  this  shape 
were  detonated. 

The  use  of  nozzle  blockage  to  modify  cloud  shape  was  abandoned 
after  these  disappointing  results.  Other  techniques  will  have  to  be 
tried  if  hemispherical  clouds  having  uniform  fuel  distribution  are 
to  be  attained.  For  example,  experiments  might  include  nozzle 
arrays  with  some  downward-directed  jets,  or  arrays  which  have 
variable  nozzle  diameters  or  jet  velocities. 
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Heptane  cloud  formed  with  dimple-type  nozzle  blockage 
pattern,  1400-nozzle  head  with  20  percent  blockage, 
test  CF-15.  The  distance  between  the  diamond-shaped 
markers  is  9.2  m  (30  ft) . 


Figure  44.  Cloud  formed  with  dimple-type  nozzle  blockage  pattern, 

1400-nozzle  head  with  10  percent  blockage,  FAE  test  CD-11, 
78.5  percent  heptane  with  21.5  percent  propyl  nitrate  by 
volume.  The  distance  between  the  diamond-shaped  markers 
is  9.2  m  (30  ft ) . 
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Figure  47.  Heptane  cloud  formed  with  ring-type  nozzle  blockage  pattern, 
1400-nozzle  head  with  10  percent  blockage,  test  CF-17.  The 
distance  between  the  diamond-shaped  markers  is  9.2  m  (30  ft). 
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Since  well-shaped  aerosol  clouds  were  not  successfully 
produced  in  these  experiments,  the  issue  of  cloud  shape  effects  on 
secondary  wave  production  in  FAE  airblasts  remains  unsettled.  For 
purposes  of  isolating  the  effects  of  cloud  shape  in  a  controlled 
manner,  it  is  suggested  that  gas-phase  fuel-air  mixtures  be 
detonated,  with  the  reactants  initially  contained  in  thin  plastic 
envelopes  (balloons)  of  various  geometric  shapes. 


IV.  PRESSURE  INSTRUMENTATION  REFINEMENTS 


The  earliest  FAE/nuclear  airblast  simulation  feasibility 
experiments  were  conducted  on  an  unimproved  test  site  area 
(References  1  and  2).  The  site  was  reasonably  level  but  it  had  not 
been  fully  graded.  It  was  located  in  a  valley  near  the  outlet  of  a 
small  canyon,  with  shallow  hills  on  two  opposite  sides. 

Pressure  transducers  in  use  at  the  time  were  Sensometrics 
Model  SP-91.  These  strain  gage  transducers  were  mounted  on  1/3  m  (1 
ft)  square  steel  or  plywood  plates  and  installed  directly  in  the 
ground  with  the  upper  surface  at  ground  level.  Output  from  the 
transducers  was  recorded  on  a  Honeywell  Model  1858  oscillograph 
having  a  frequency  response  to  5  kHz.  In  most  of  the  tests, 
however,  the  pressure  transducer  signals  were  filtered  with  a 
first-order  lowpass  circuit  which  limited  the  effective  system 
response  to  about  1-2  kHz.  Filtering  was  incorporated  to  eliminate 
intractable  higher  frequency  noise,  and  although  this  measure  also 
rounded  shock  wave  peaks  and  obscured  some  waveform  details,  those 
details  were  not  considered  essential  to  the  objectives  of  the  early 
efforts . 

Examples  which  characterize  the  data  taken  from  early  FAE 
tests  are  given  in  Figures  48  and  49  (Reference  19).  These  samples 
are  taken  from  the  same  experiment.  The  two  transducers  were 
located  at  approximately  the  same  radial  distance,  but  they  were  90 
degrees  apart.  A  number  of  impulsive  pressure  excursions  are  evi¬ 
dent  on  both  waveforms,  and  in  addition,  what  appears  to  be  a  sig¬ 
nificant  secondary  wave  is  present  on  Figure  48.  The  substantial 
rotational  asymmetry  between  pairs  of  pressure  traces  such  as  these 
was  partly  responsible  for  uncertainty  as  to  whether  the  waveform 
aoerrations  were  a  result  of  peculiarities  inherent  to  fuel -air 
explosions,  or  whether  these  disturbances  could  be  attributed  to 
cnaracteristics  of  the  instrumentation  layout  and  installation. 


It  was  accordingly  undertaken  to  improve  the  test  site,  with 
careful  attention  given  to  the  transducer  mountings  and  to  their 
locations  with  respect  to  adjacent  hills.  The  upgraded  facility  is 
described  in  detail  in  Reference  2.  At  the  same  time,  piezoelectric 
pressure  transducers  (PCS  Model  102A12)  were  acquired  to  replace  the 
strain  gage  pressure  transducers.  In  subsequent  FAE  testing,  the 
piezoelectric  transducers  were  used  in  conjunction  with  an  EMI  Model 
SE7000A  14-channel  tape  deck.  At  a  recording  velocity  of  305  cm/s 
(120  in/s ) ,  this  tape  recorder  has  a  frequency  response  to  80  kHz. 
It  was  the  limiting  component  in  the  new  data  acquisition  chain. 
Filtering  was  not  used  on  the  input  to  this  device. 

During  the  first  tests  of  the  improved  experimental  facility 
and  instrumentation*  (FAE  tests  CD-I  to  CD-5),  the  pressure  records 
from  transducers  close-in  to  the  blast  source  exhibited  considerable 
ringing.  A  worst-case  example  is  given  on  Figure  50.  The  ringing 
was  not  a  simple  sinusoid;  it  could  be  ameliorated  but  not 
eliminated  by  a  data  processing  procedure  involving  narrowband 
filtering  of  the  waveform  Fourier  transform.  It  was  felt  that  this 
ringing  problem  compromised  the  virtues  of  the  instrumentation 
improvements,  and  that  further  FAE  experiments  should  be  postponed 
until  the  problem  was  resolved. 

Piezoe’ ectr ic  pressure  transducers  respond  to  acceleration  as 
well  as  to  pressure.  Generally,  the  extent  of  sensitivity  to 
acceleration  varies  with  the  direction  of  movement.  Most 
commercially  available  transducer  designs  include  components  which 
partially  compensate  at  least  for  axial  acceleration  sensitivity. 
(A  comparison  of  some  widely  used  airblast  pressure  transducers 
appears  on  Table  2.)  It  was  believed  that  physical  oscillations  of 
the  pressure  transducers  were  responsible  for  the  pressure  waveform 
ringing. 

A  number  of  steps  were  taken  to  refine  the  pressure 
instrumentation  so  as  to  reduce  ringing  and  improve  data 
confidence.  These  steps  included  modifications  to  the  transducer 


*  See  Appendix  A. 
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Figure  50.  Pressure  record  obtained  from  FAE  test  CD-5,  0.265  m3  (7  qal . )  propylene  oxide 
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TABLE  2 

PRESSURE  TRANSDUCER  PERFORMANCE  SPECIFICATIONS 
(DATA  DERIVED  FROM  MANUFACTURERS  BROCHURES) 


Transducer 

PCB 

Endevco 

Kullte 

Klstler 

K1 stler 

Type 

102K58 

Piezoelectric 

8510 

Plezoresl stive 

XTE-190 

Plezoresistlve 

603B 

Piezoelectric 

Built  In  amp. 

Charge  Mode 

Built  In  amp. 

Range 

500  psi 

15  psi -500  psi 

25  psi -500  psi 

30  psi -3000  psi 

100-500  psi 

Resonant 

Frequency 

(KHz) 

500 

100-500 

100-350 

>400 

130 

Rise  time 
(usee) 

1.0 

— 

— 

1.0 

3.0 

Tine  Constant 
(sec) 

>5 

~  30 

— 

25-110 

Acceleration 

Sensitivity 

(psl/g) 

(Axial ) 

0.002 

0.0004 

0.0001-0004 

<  0.0015 

0.015 

Flash 

Sensitivity 

(ps1/*F) 

0.003-0.01 

— 

Sensitivity 
(mv/psl ) 

+20 

+20-+0.63 

+20-+-.15 

+C.35/pc/psi 

-50-10 

LI nearl ty 
(+  percent 
FTS. ) 

1.0 

0.5 

1.0 

1.0 

mountings,  construction  of  a  device  for  in  situ  pressure 

calibration,  and  testing  of  transducer-face  thermal  shields. 

1.  TRANSDUCER  MOUNTING 

A  cross-sectional  sketch  of  the  original  transducer  mounting 
arrangement  is  given  in  Figure  51.  The  sketch  is  drawn  to  scale. 
Transducers  were  installed  directly  into  the  steel  mounting  plates, 
wnich  have  dimensions  of  15.2  cm  x  61  cm  x  1.9  cm  (6  in.  x  24  in.  x 
3/4  in.).  The  cover  plates  were  secured  to  steel  angle-iron 
supports  with  three  recessed  socket  head  cap  screws  on  ^  f  the 
two  long  sides.  This  arrangement  was  used  in  FAE  tests  to  CD-5. 

In  a  first  attempt  to  reduce  ringing,  the  mou^ti  -,s  were 
modified  by  installing  the  transducers  in  simple  plastic  inserts 
whicn  were  epoxied  into  the  steel  plates.  Delrin  inserts  were  tried 
on  FAE  test  CD-6,  and  teflon  inserts  were  tried  on  the  next  and 
otherwise  identical  test,  CD-8.  The  transducer  ringing  was  not 
suppressed  in  either  of  these  two  tests. 

A  special  FAE  test  ( COP— 1 )  was  then  conducted  in  which  four 
additional  transducer  mounting  variations  were  examined  for  their 
effectiveness  in  reducing  ringing.  The  arrangements  tried  were  as 
fol lows: 

a.  A  transducer  was  installed  in  a  threaded  washer 
which  was  potted  into  Biwax  Corporation  Type 
601  encapsulant,  as  shown  schematically  in 
Figure  52(a).  This  encapsulant  has  been  found 
in  other  applications  to  very  effectively  damp 
mechanical  vibrations. 

b.  A  transducer  was  potted  into  Biwax  601  as  in 
(a),  but  the  mounting  plate  was  laminated.  The 
laminations  comprised  a  top  layer  of  7.94  mm 
(5/16  in.)  steel  and  a  bottom  layer  of  6.35  mm 


strument 


ss-sectional  view  of  instrumentation  trough 


(1/4  in.)  steel.  A  4.76  mm  (3/16  in.)  thick 
layer  of  a  viscoelastic  compound  (E.A.R. 

Corporation  type  C2003)  was  glued  between  the 

steel  plates  with  Bostik  adhesive.  This 
arrangement  is  sketched  in  Figure  52(b).  The 

overall  thickness  of  this  lamination  was  the 
same  as  the  thickness  of  the  other  plates. 

c.  A  mounting  plate  was  made  in  which  lamination, 

as  in  (b)  above,  was  used,  but  without  the 

Biwax  suspension.  The  transducer  was  installed 
directly  in  the  top  steel  plate  of  this 
laminate. 

d.  One  of  the  mounting  plates  was  rigidly  pinned 
to  the  concrete  floor  of  the  instrumental  ion 
trough  in  which  it  was  installed.  The  plate 
was  securely  bolted  to  five,  5  cm  (2  in.)  dia. 
steel  reinforcing  rods  that  had  been  grouted 
into  the  instrumentation  trough  with  Burkestone 
hydraulic  concrete.  The  plate  itself  was  made 
of  plain  steel  with  a  Oelrin  insert  supporting 
the  transducer. 


In  addition  to  these  four  mounting  arrangements,  one 
transducer  was  directly  installed  on  a  plain  steel  plate.  On  this 
same  plate,  a  dummy  transducer  was  installed  in  a  blind  hole.  All 
of  the  plates  were  located  on  the  test  pad  within  13.9  m  (45.6  ft) 
of  the  explosion  center.  (Transducer  ringing  had  not  previously 
been  significant  beyond  this  range.) 

When  test  CDP-1  was  conducted,  none  of  the  four  mountinq 
techniques  proved  to  be  effective  in  significantly  reducing 
ringing.  The  mounting  arrangement  in  which  both  Biwax  suspension 
and  plate  lamination  were  used  yielded  the  best  results,  but  ringing 
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was  not  eliminated.  The  results  from  the  plain  steel  plate  (on 
which  both  active  and  dummy  gages  were  mounted)  appear  on  Figure 
53.  The  major  frequency  component  in  the  pre-shock  ringing  of  the 
active  transducer  was  about  400  Hz.  The  major  frequency  component 
in  tne  dummy  gage  signal  (which  has  been  scaled  to  pressure  units) 
is  about  1500  Hz. 

In  a  second  FAE  test  of  this  type  (designated  CDP-2),  pairs  of 
active  transducers  were  mounted  directly  on  each  of  three  plain 
steel  plates.  Each  pair  consisted  of  a  PCB  piezoelectric  transducer 
and  an  Endevco  Model  8510  piezoresistive  strain  gage  transducer. 
The  resonant  frequencies  of  these  devices  are  500  kHz  and  240  kHz, 
respecti vely. 

The  objective  of  CDP-2  was  to  evaluate  the  relative 
performance  of  these  two  transducers  when  mounted  in  the  simplest 
fashion.  The  Endevco  transducer  was  chosen  for  comparison  with  the 
PCB  because  it  has  equivalent  high  frequency  response  capability, 
but  is  of  a  different  generic  type,  with  much  lower  acceleration 
sensitivity.  Figure  54  gives  the  results  of  drop  bar  tests  of  these 
two  transducers,  using  the  apparatus  described  in  References  20  and 
21. 

Figure  55  shows  the  results  obtained  from  one  of  the  three 
transducer  pairs  during  test  CDP-2.  The  two  pressure  records  in 
this  pair  have  been  superimposed  in  absolute  time  with  respect  to 
detonation  initiation.  The  overshoot  that  can  be  seen  at  the 
leading  shock  appeared  only  on  the  Endevco,  a  character istic  of  that 
transducer  which  was  also  evident  in  shock  tube  tests.  The  Endevco 
record  appears  somewhat  more  finely  detailed  than  that  obtained  from 
the  PCB,  possibly  due  to  the  difference  in  sensing  element  size. 
The  sensing  element  diameter  of  the  PCB  is  approximately  5.8  mm 
(0.23  in.),  while  that  of  the  Endevco  is  2.5  mm  (0.10  in.). 

On  Figure  55  there  appears  to  have  been  some  negative 
displacement  of  the  average  PCB  signal  level.  This  is  presently 
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Active  dnd  dummy  pressure  gage  records  from  FAE  test  CDP-1,  0.0265  m3  (7  gal.) 
propylene  oxide,  R  =  9.93  m  (32.6  ft.).  The  dummy  gage  output  has  been  converted 
to  virtual  pressure  units  for  comparative  purposes. 
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Figure  54.  Results  of  drop  bar  acceleration  sensitivity  tests  of  PCB  and  Endevco  transducers 
Peak  acceleration  was  ^980  m/s2  (100  g)  in  both  cases.  Gage  output  was  converted 
to  virtual  pressure  units  for  comparitive  purposes. 


ENDEVCO 


Pressure  records  obtained  from  FAE  test  C 
Both  transducers  mounted  at  same  location 


oelieved  to  be  due  to  inadequate  thermal  protection  of  the  sensing 
element  at  the  time  of  this  test.  The  sensing  element  of  the 
Endevco  was  protected  with  a  graphite/silicon-grease  compound  (on 
the  recommendation  of  the  manufacturer).  However,  the  coating  used 
on  the  PCB  was  plain  silicon  grease,  which  was  later  shown  to  be 
ineffective  as  a  thermal  shield. 

Noticible  pre-shock  ringing  appeared  only  on  the  PCB  pressure 
record.  There  was  also  a  rather  significant  quantitative  dis¬ 
crepancy  between  the  two  records  during  the  first  half-millisecond 
following  shock  passage.  However,  apart  from  the  apparent  negative 
drift  in  the  PCB  signal,  the  two  transducers  tracked  very  closely 
throughout  the  remaining  positive  phase. 

Considering  the  magnitude  of  difference  in  acceleration 
sensitivity  between  the  PCB  and  Endevco  transducers  that  were  used 
in  test  CDP-2,  it  was  concluded  from  this  that  the  post-shock 
oscillations  on  the  pressure  waveforms  in  Figure  55  must  have  been 
largely  of  gasdynamic  origin.  It  was  noted  that  the  lower  amplitude 
pressure  oscillations  were  most  prevalent  in  those  records  from 
transducer  locations  that  are  presumably  engulfed  by  the  expanding 
fuel-air  cloud  combustion  products.  It  was  felt  that  some  of  the 
oscillations  at  these  locations  may  have  been  due  to  combustion 
product  turbulence.  In  any  event,  the  oscillations  did  not  appear 
to  be  entirely  the  result  of  transducer  acceleration  sensitivity. 

Since  transducer  ringing  was  not  alleviated  by  simple 
mod  1  f icat ions  in  the  mounting  arrangements,  the  decision  was  made  to 
circumvent  tne  proolem  in  future  work  by  employing  piezores i s t i ve 
t/pe  transducers  at  those  locations  that  are  close  to  the  explosion 
source .  Ihe  tradeoffs  which  may  have  to  be  made  include  the 

o« er snoot  pn. olein  with  some  piezoresistive  transducers,  as  well  as 
posjihie  difficulties  with  durability  and  serviceability. 

•* ter»Mt  i  vc  1 1 ,  toe  use  of  piezoelectric  transducers  might  be 
•■on tinucu,  out  individual  units  should  then  be  carefully  selected 
r  ,r  tneir-  low  acceleration  sensitivity.  In  drop  tests  of  those  PCB 


transducers  that  were  used  in  the  presently  reported  FAE 
experiments,  acceleration  sensitivities  in  terms  of  virtual  pressure 
varied  over  more  than  an  order  of  magnitude,  from  a  minimum  of  1.4 
Pa-s  /m  (0.002  psi/g)  to  a  maximum  of  21  Pa-s  / m  (0.03  psi/g). 

2.  IN-SITU  PULSE  PRESSURE  TRANSDUCER  CALIBRATOR 

Most  commercial  pressure  transducers  are  subject  to  changes  in 
calibration  due  to  weather  conditions,  aging  of  parts,  connector 
cleanliness,  and  cable  properties.  Hysteresis  effects  may  alter  the 
calibration  during  loading.  In  order  to  establish  confidence  in  the 
accuracy  of  pressure  measurements,  in  situ  calibration  is  essential, 
preferably  before  and  after  each  experiment.  For  purposes  of 
calibrating  piezoelectric  type  transducers,  a  device  is  needed  that 
will  apply  pressure  to  the  sensing  element  in  a  time  period  that  is 
short  compared  to  the  decay  time  constant  of  the  transducer. 

A  pulse  calibrator  of  this  type  was  designed  and  built  for 
in-situ  use  in  conjunction  with  the  FAE  experiments.  A  photograph 
of  the  device  appears  in  Figure  56.  In  this  figure,  a  static 
pressure  transducer  is  shown  installed  in  the  steel  cover  plate  of 
one  of  the  instrumentation  troughs.  The  pulse  calibrator  head  is 
resting  on  the  cover  to  the  right  of  the  transducer.  In  use,  the 
head  would  be  centered  over  the  transducer  by  means  of  scribed 
alignment  markings  on  the  cover  plate. 

The  calibrator  head  is  a  modified  D.C.  electromagnet  (Magnet 
Sales  and  Manufacturing  Model  EM-R6B)  having  a  holding  force  of  2670 
N  (6U0  lb).  An  0-ring  groove  was  machined  into  the  face  of  this 
magnet.  A  small  circular  recess,  12.7  mm  (1/2  in.)  dia.  x  1.59  mm 
(1/16  in.)  deep,  was  cut  into  the  center  of  the  magnet  face  to 
prevent  mechanical  contact  with  the  transducer,  and  a  3.18  mm  (1/8 
in.)  diameter  hole  was  drilled  through  the  axis  of  the  magnet.  An 
electric  solenoid-actuated  vulve  (Skinner  type  X54)  was  mounted  as 
closely  as  possible  to  the  top  surface  of  the  magnet.  The  solenoid 
valve,  which  normally  operates  on  24  VDC,  has  an  effective  orifice 


Figure  56.  In-situ  pressure  transducer  pulse  calibrator. 
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diameter  of  4.76  mm  (3/16  in.),  and  a  nominal  opening  time  of  16 
ms.  For  this  application,  faster  opening  was  obtained  by 
discharging  a  3000uF  capacitor  at  250V  through  the  solenoid  with  a 
momentary  switch.  The  current  decays  to  a  steady-state  level 
consistent  with  continuous  solenoid  operation  through  the  use  of  a 
resistor  in  series  with  the  250V  battery.  The  decay  time  constant 
is  on  the  order  of  200  ms. 

The  pressure  source  for  calibration  is  nitrogen  gas,  which  is 
stored  in  a  0.0038  m^  (1  gal)  bottle  built  into  the  portable  gage 
case  shown  behind  tne  instrumentation  trough  in  Figure  56.  This 
bottle  is  initially  charged  from  a  gas  cylinder  to  the  highest 
pressure  needed  in  the  calibrations.  Nitrogen  is  then  exhausted 
from  the  bottle  as  needed  for  lower  pressure  calibrations.  The 
bottle  is  connected  to  the  calibrator  head  with  6.4  mm  (1/4  in.) 
I.D.  plastic  tubing.  The  very  small  volume  expansion  of  the 
calibration  gas  that  accompanies  opening  of  the  solenoid  valve 
assures  that  tne  drop  in  source  pressure  will  be  negligible. 

The  source  pressure  is  measured  by  one  of  four  Ashcroft  Model 
i279A  Bourdon  tube  pressure  gages.  These  gages  have  a  full  scale 
accuracy  of  +  0.5  percent.  Their  accuracy  at  part  scale  is 

proportionately  lower.  However,  the  ranges  of  the  four  gages  used 
with  the  calibrator  were  chosen  with  sufficient  overlap  so  that  the 
maximum  error  of  any  reading  over  the  entire  pressure  range  of 
interest  would  not  exceed  approximately  +  1.5  percent. 

In  order  to  calibrate  stagnation  pressure  stings  in  situ,  the 
adaptor  rack  shown  in  Figure  57  was  built.  The  sting  is  inserted 
into  a  retainer  that  was  machined  to  fit  the  contour  of  the  sting 
tip.  A  Conax  packing  gland  type  feedthrough  connector  securely 
nolds  the  sting  in  position,  with  its  active  surface  flush  with  the 
steel  inset  on  the  adaptor  rack.  The  sting  is  shown  installed  in 
the  adaptor  in  Figure  57. 


9 


114 


i 


When  the  holding  magnet  is  initially  energized,  a  small 
transient  signal  is  generated  by  the  transducer.  This  signal  is 
believed  to  be  due  to  three  mechanical  effects:  a  slight 

displacement  of  internal  components  in  the  transducer  by  the 
magnetic  force,  impact  of  the  electromagnet  face  with  the  trough 
cover,  and  a  slight  increase  of  pressure  in  the  volume  above  the 
transducer  due  to  the  decrease  in  this  volume.  The  transient  signal 
decays  quickly  and  does  not  interfere  with  the  subsequent 
calioration  process. 

An  example  of  the  performance  of  the  pulse  calibrator  is  given 
in  Figure  58.  In  this  example,  the  supply  pressure  was  0.69  MPa 
(100  ps i ) .  The  pressure  risetime  was  about  7  ms.  Oscillations  with 
a  mean  frequency  of  about  80  Hz  and  having  a  half-wave  amplitude  of 
4  percent  of  the  average  signal  level  continued  for  approximately 
Idu  ms.  The  present  overall  calibration  accuracy  is  therefore 
regarded  as  better  tnan  4  percent.  It  may  be  possible  to  improve 
tms  accuracy  somewhat  by  reducing  the  initial  oscillations,  through 
tne  use  of  a  solenoid  valve  having  a  somewhat  smaller  effective 
orifice  size  (Reference  22). 

The  in-situ  calibrator  is  useful  for  testing  blast  type 
pressure  transducers  immediately  before  and  after  FAE  experiments. 
However,  for  assessment  of  response  to  shock  waves,  the  transducers 
are  periodically  tested  in  a  shock  tube.  An  in-situ  shock  tube 
calibrator  of  the  type  described  in  Reference  23  could  also  be  used 
for  this  latter  purpose. 

3.  THERMAL  SHIELDING 

The  piezoelectr ic  pressure  transducers  used  in  the  FAE  experi¬ 
ments  were  found  to  be  very  sensitive  to  transient  heating.  The 
magnitude  of  error  introduced  into  the  positive-phase  FAE  blastwave 
pressure  records  by  this  effect  was  small  but  noticible  at  the 
1/4-ton  (nominal  nuclear  equivalent)  scale.  At  that  scale,  most  of 
the  positive-phase  durations,  measured  at  ranges  of  interest,  are 


less  tnan  50  ms;  those  at  positions  having  the  greatest  thermal 
loading  are  less  than  about  20  ms. 

Figure  59(a)  shows  a  200  ms  static  pressure  record  from  an  FAE 
test,  at  a  location  that  was  just  outside  the  initial  cloud  radius. 
The  active  surface  of  the  transducer  was  in  this  case  coated  with  a 
thin  film  of  plain  silicone  grease.  On  other  FAE  tests,  a  single 
layer  of  0.178  mn  (0.007  in.)  thick  black  plastic  electrical  tape 
had  previously  been  used  as  a  thermal  and  ablative  shield,  with 
similar  results.  The  dramatic  improvement  in  thermal  shielding 
ootained  by  adding  zinc  oxide  to  the  silicone  grease  coating,  and 
used  in  conjunction  with  a  single  layer  of  black  electrical  tape,  is 
illustrated  oy  the  example  in  Figure  59(b).  The  two  traces  in  this 
figure  are  from  different  tests  and  taken  with  different  transducers 
(of  the  same  manufacture] .  However,  the  test  conditions  were  very 
similar,  and  the  transducers  were  located  in  the  same  position  on 
tne  test  pad. 

Several  types  of  easily  applied  thermal  shields  were  systemat¬ 
ically  tested  following  the  procedure  described  in  Reference  24.  A 
shielded  test  transducer  and  an  unshielded  control  transducer  were 
flush  mounted  on  a  metal  plate.  At  a  distance  of  70  rrm  (2.75  in.) 
along  the  normal  to  the  surface  of  this  plate,  a  No.  22  photographic 
flash  bulb  was  fired.  This  produces  radiative  rather  than  convec¬ 
tive  tnermal  loading,  but  it  is  felt  that  the  relative  effects  of 
various  forms  of  shielding  would  be  similar  in  either  case. 

The  output  from  the  test  transducer  using  different  shield 
types  is  shown  on  Figure  60.  Tne  most  effective  shield  tested 
consisted  of  a  single  layer  of  olack  electrical  tape,  coated  with  a 
commercial  zinc-oxide/silicone  grease  heat  sink  compound  (Chem- 
tronics  SL-lJ.  Three  layers  of  black  tape  (not  shown  on  Figure  60) 
were  almost  as  effective,  but  there  was  concern  that  this  amount  of 
tape  might  interfere  with  tne  transducer  response. 
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(a) 


(b) 

Figure  59.  Pressure  records  obtained  from  FAE  tests,  0.0265  m3  (7  gal.) 
propylene  oxide,  R  =  5.92  m  (19.4  ft.),  (a)  Test  CDP-2, 

transducer  coated  with  silicon  grease,  (b)  Test  CD-14,  trans¬ 
ducer  covered  with  black  plastic  electrical  tape  and  coated 
with  sil iconegrease/zinc  oxide  compound. 
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At  present,  a  single  layer  of  electrical  tape  coated  with  heat 
sink  compound  is  used  with  all  piezoelectric  transducers  for  FAE 
airblast  measurements.  The  active  face  of  the  transducer  is  cleaned 
with  acetone  and  the  tape  is  carefully  applied  to  avoid  trapped  air 
pockets.  The  transducer  is  then  installed  in  a  mounting  plate 
(Figure  51).  A  jig  is  used  so  that  when  the  transducer  is  securely 
fastened  in  place,  the  upper  surface  of  the  tape  is  0.635  mm  (0.025 
in.)  below  the  plate  surface.  The  cavity  above  the  transducer  is 
then  filled  with  the  heat  sink  compound  until  the  compound  is  flush 
with  the  top  of  the  plate. 


V.  CONCLUDING  REMARKS 


Small-scale  fuel-air  clouds  have  been  impulsively  formed  from 
a  central  source  and  detonated,  but  these  clouds  are  not  hemi¬ 
spherical,  the  positive-phase  pressure  decay  is  not  always  mono¬ 
tonic,  and  the  pressure  profiles  frequently  exhibit  secondary  wave 
systems.  Large-scale  impulsive  jets  from  single  nozzles  have  been 
produced,  but  the  reach  needed  for  a  1  KT  simulator  has  not  yet  been 
attained.  The  data  needed  to  realistically  assess  technical  feas¬ 
ibility  of  a  full-scale  simulator  installation  have  not  been 
acquired.  To  complete  the  analysis  of  feasibility  for  a  1  KT 
FAE/nuclear  airblast  simulator,  a  number  of  additional  investi¬ 
gations  are  needed. 

a.  Physical  Feasibility 

An  important  remaining  question  is  whether  secondary  shocks 
and  other  waveform  aberrations  can  be  eliminated  from  the  FAE 
airblast.  If  these  are  the  result  of  cloud  shape  distortions  then 
the  question  is  more  of  a  technological  issue  than  a  basic  question 
of  the  FAE  physics.  Significant,  impulse-carrying  secondary  waves 
were  not  clearly  evident  in  the  balloon  tests  conducted  as  part  of 
Operation  Distant  Plain  (References  8  and  9).  In  those  tests, 
vapor-phase  fuel  was  mixed  with  oxygen  and  the  mixture  confined  in  a 
balloon  envelope  prior  to  detonation.  This  indicates  that,  under 
conditions  of  accurate  cloud  shape,  reasonable  mixture  uniformity, 
and  point  symmetric  initiation,  significant  secondary  wave  systems 
may  not  form  or  may  be  essentially  inconsequential.  The  simplified 
S3  FAE  airblast  calculation  (Reference  2)  also  showed  secondary 
waves  only  in  the  negative  phase  under  these  conditions.  However, 
that  calculation  assumed  that  the  combustion  products  were 
chemically  frozen.  To  investigate  the  possibility  that  second 
shocks  might  result  from  rapid  chemical  composition  changes  induced 
by  cloud  expansion,  this  calculation  could  be  extended  to  include 
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■  euji  equilibrium  and  variable  specific  heats.  The  effects  on 
jc^onoary  wave  formation  of  cloud  shape  asymmetry,  initiator 
location  asymmetry,  an  extended  detonation  reaction  zone,  and 
f  je 1 / a  i  r  mixture  inhomogeneities  could  also  be  treated 
ana  1 yt ica 1  ly .  Additionally,  small-scale  FAE  balloon  tests  could  be 
conducted  to  extend  the  Distant  Plain  data,  with  the  objective  of 
examining  the  existence  of  secondary  waves  under  optimal 
conditions.  Many  of  the  problems  with  large-scale  balloon  testing 
are  more  manageable  at  smaller  scales. 

An  expanded  FAE  data  base  is  also  needed  to  enable  final 
judgements  to  be  made  as  to  the  quality  of  simulation,  repeatability 
and  symmetry  of  the  airblast  that  is  generated  when  an  impulsively 
formed  fuel-air  cloud  is  detonated.  The  data  base  should  be  updated 
as  improvements  in  cloud  shaping  and  fuel  distribution  uniformity 
are  accomplished. 

b.  Technological  Feasibility 

Cloud  formation  remains  one  of  the  most  serious  technological 
problems  that  could  impede  the  development  of  a  reusable  FAE/nuclear 
airblast  simulator.  The  airblast  from  a  non-hemi spherical  FAE  is 
not  likely  to  be  point  symmetric  near  the  exploded  cloud.  It  also 
appears  possible  that  cloud  misshapings  may  contribute  to  secondary 
wave  systems.  Further  efforts  in  the  area  of  cloud  formation 
technology  therefore  seem  essential  to  the  advancement  of  the 
overall  feasibility  assessment. 

c.  Engineering  Feasibility 

Most  of  the  previous  experimental  work  dealing  with  single, 
impulsive  jets  has  been  of  a  survey  nature.  Systematic  experiments 
at  small  and  intermediate  scales  are  now  needed  in  order  to  more 
accurately  anticipate  the  minimum  requirements  of  a  large-scale 
dispensing  system.  These  tests  should  include  the  following 
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variables:  nozzle  diameter,  dispensing  pressure,  volume  of  fuel 
dispensed,  constancy  of  dispensing  pressure,  and  the  effects  of  fuel 
properties  and  additives.  In  addition,  optional  fuel  pressurization 
schemes  should  be  examined  at  the  smallest  practical  scale. 
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D  USING  UPGRADED  FACILITIES 


Propyl  nitrate  7.  Tests  C01-C08  conducted 

under  separate  contract 
(Ref.  2) 


APPENDIX  B. 


LISTING  OF  TESTS  CONDUCTED  TO  EXAMINE  CLOUD  SHAPE  MODIFICATIONS 

WITH  NOZZLE  BLOCKAGE  U) 


Test(5) 

Date 

Fuel(3) 

Dispensing 

Pressure 

(psi) 

Blockage^ 

Type 

Percent 

CF16 

10-2-80 

Heptane 

244 

Dimple 

2 

CF20 

10-13-80 

Heptane 

232 

Dimple 

10 

CD11 

10-15-80 

Hept/PN^ 

218 

Dimple 

10 

CD12 

10-22-80 

Hept/PN^ 

220 

Dimple 

15 

CF15 

9-25-80 

Heptane 

270 

Dimple 

20 

CDIO 

9-23-80 

Hept/PN^ 

250 

Dimple 

25 

CF18 

10-8-80 

Heptane 

204 

Ring 

5.5 

CF19 

10-10-80 

Heptane 

240 

Ri  ng 

7.5 

CF17 

10-7-80 

Heptane 

243 

Ring 

10 

NOTES: 

1.  Nozzle  head  has  1400  0.116  in. 

(4)  Heptane  and  Propyl  nitrate 

dia.  holes  (w/o  blockage) 

(See  Appendix  A) 

• 

2.  See  text 

(5)  CF  Designates  test  for  cloud 

formation  only  (no  detonation); 

3.  7  gallons  dispensed  in 

CD  designates  cloud  detonation 

all  tests 

test 

« 


DEPARTMENT  GF  NAVY  4  Conti nued ) 

*  r'd  te.;  it  Systems  Pro1  net  Office 
ATTN:'  NbP-y/L 
A  TIN:  N  S  P  -  A  3 

or  PARENT  OP  THE  AIR  F  ORCL 


ATTN:  I  NT 

Air  ’01\>J  GeOph  ys  ICS  Laboratory 
ATTN:  Lwh,  h.  Ossirg 

Atr  Institute  of  Technology 

ATTN:  Library 

inters 

A.  •  * .  i  r ,  e  S  .  s  terrs  C orana  nd 

A  i  *  f.r.t-  rteapivs  Laboratory 
A T TN  :  NTT,  M .  PI  jn.ot ide n 
ATTN  :  N  T  E  D ,  i  .  Seusy 
ATTN:  NT t  S 
^TTN:  SuL 
A 'TV  NTT  0 

'*  i  r  .  n  ’  v  e»*s  i  t.  >  Library 
A ’  T N :  A.,l  -LSt 

a  I  i  n  :  i  t.  M  l  s  i  I  e  Of  f  i  i  <.•  DAA 
ATTN;  Pp 

e . T N :  [ NOT  ,  0.  Ga Oe 

.  L ,  ATTN:  ENSN 

Jt'Pj  t  .  'Jl  1 1‘*'  U*  Staff 

'■e^eir s h,  Development,  N  Acq 
ATTN:  AFRDQI 

■»*r.-ito;i<  Air  Comma n,c! 

ATTN:  DO TP 

ATTN:  NR  I  ST  INFO  Library 
l  PA.kTM:  NT_  OF  ENERGY 

•  •  r  *  ■  vn  T  1  *  L  ner  ;v 

•It-.  e’roMt*  Opera  *  i  Vo.  Of  f  ice 


’TN:  0.  Richmond 

'*  v  *  ■  t  ! r ,  Apcluat ion 


;  •ear  t. on-.  Office 

"N:  i  oo  Con  for  N  chnkal  Library 

. :  AN Vj;  NT  AGENCY 

I1  Tell i qerce  Agency 
*N:  US*R  NT D 


A  "  N  :  ■ ;  *,  Do  t.  j;,;e r.  *.  s  0  f  *  i » .  •  ■  r 

1  Ak  T  vt  N  T  J  I.NLPijY  CONTRAC 'OR S 

:  india  Na*i  cn.il  Labor.*  *  or  lr** 

AT  'N  :  l  i  b  S  Sei.  1  a  -s  Di  v 


DEPARTMENT  Of  ENERGY  CONTRACTOR^  _{ Continued j 

University  of  California 
Lawrence  Livermore  National  Lab 

ATTN:  Technical  Info  Dept.  Library 

1  os  Alamos  National  Laboratory 
ATTN;  Reports  Library 
ATTN:  J.  Hopkins 

Sandia  National  Laboratories 
ATTN:  Tech  Lib  3141 

DEPARTMENT  0  M)EF IN  IS  E  __C_0_N  T  R  ACT  0  R  S 

Aerospace  Corp 

ATTN:  Library  Acquisition  Ml  199 

Applied  Research  Associates,  Inc 
ATTN:  N.  Higgins 
ATTN :  J.  Bratton 

Applied  Research  Associates,  Inc 
ATTN;  J.  Shinn 

Applied  Research  Associates,  Inc 
ATTN:  D.  Piepenburg 

Applied  Research  Associates,  Inc 
ATTN:  R.  Frank 

ROM  Corp 

ATTN:  Corporate  Lib 
ATTN:  T.  Neighbors 

Boeing  Co 

ATTN:  Aerospace  Library 

California  Research  &  Technology,  Inc 
ATTN:  K.  Kreyenhagen 
ATTN:  M.  Rosenblatt 

California  Research  A  Technology,  Inc 
ATTN:  F.  Sauer 
ATTN:  D.  Orpha ] 

Cushing  Associates 

ATTN:  V.  Cushing 

F.G&G  Wash  Analytical  Svcs  Ctr,  Inc 
ATTN:  Library 

E lectro-Mech  Systems,  Inc 
ATTN:  R.  Shunk 

EIoctrn-Mech  Systems,  Inc 
ATTN:  L.  Piper 

GEO  Centers,  Inc 

ATTN:  E.  Mar  rani 

H-Toch  Labs,  Inc 

A1TN:  B.  Ha»'tenbau!’- 

Ill  Resear  h  InstHjfe 

ATTN:  Documents  Library 

institute  for  Defense  Analyst.-- 
ATTN:  Cl  ass  i  *  i**d  Library 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

,'AYCOR 

ATTN: 

Library 

Raman  Sciences  Corp 

ATTN: 

Library 

samar,  Sciences  Corp 

ATTN: 

E.  Conrad 

kasan  Tempo 

ATTN: 

DAS! AC 

Karan  Tempo 

ATTN: 

DAS  I  AC 

Lockheed  Missiles  X  Space  Co,  In c 

ATTN: 

T.  Geers,  Orgn  52-33 

ATTN: 

Technical  Library 

Locuneed  Missiles  X  Space  Co,  Inc 

ATTN: 

S.  Taimuty,  Dept  81-74/154 

ATTN: 

TIC-Library 

Martin  Marietta  Denver  Aerospace 

ATTN: 

D-6074 ,  G.  Freyer 

witre  Corp 

ATTN: 

MS  E19Q 

National  Technical  Systems 

ATTN: 

P.  'eberman 

„n> vers i ty 

of  New  Mexico 

ATTN 

G.  Leigh 

ATTN 

N.  Baum 

i  cy  ATTN 

D.  Calhoun 

,.ni  vers )  t  v 

of  New  Mexico 

ATTN 

CERF,  N,  Baum 

Pacific-Sierra  Research  Corp 

ATTN 

H.  Brode,  Chairman  SAGE 

Pacifica  Technology 

•MIN 

G.  Kent 

in, sics  Applications,  Inc 

ATTN 

f.  Ford 

i’t  International  Co 

ATTN 

C.  Moore 

ATTN 

Technical  Library 

A Associates 

AT  IN 

:  J.  Lewis 

ATTN 

:  Technical  Information  Center 

ATTN 

D.  Simons 

ATTN 

:  C.  Lee 

ATTN 

P.  Maas 

DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 
Rand  Corp 


ATTN: 

P.  Davis 

Rand  Corp 

ATTN: 

B.  Bennett 

S-CUBED 

ATTN: 

D.  Grine 

ATTN: 

Library 

2  cy  ATTN: 

T.  Pierce 

2  cy  ATTN: 

M.  Groethe 

2  cy  ATTN: 

R.  Sedgwick 

Science  Applications,  Inc 

ATTN:  Technical  Library 

Science  Appl  ications,  Inc 
ATTN:  M.  Knasel 
ATTN:  W.  Layson 
ATTN:  J.  Cockayne 

Southwest  Research  Institute 
ATTN:  A.  Wenzel 
ATTN:  W.  Baker 

SRI  International 

ATTN:  J.  Colton 
ATTN:  G.  Abrahamson 

Structural  Mechanics  Associates,  Inc 
ATTN:  R.  Kennedy 

Teled/ne  Brown  Engineering 
ATTN:  F.  Leopard 
ATTN:  D.  Ormond 

TRW  Electronics  X  Defense  Sector 

ATTN:  Technical  Information  Center 
ATTN:  D.  Baer 
2  cy  ATTN:  N.  Lipner 

TRW  Electronics  X  Defense  Sector 
ATTN:  P.  Dai 
ATTN:  E.  Wong 

We  id  I inger  Assoc,  Consulting  Engrg 
ATTN:  M.  Baron 

Weidl inger  Assoc,  Consulting  Engrg 
ATTN:  J.  Isenherg 

Weidl inger  Assoc,  Consulting  Engrg 
ATTN:  A.  Misovec 
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